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A MATHEMATICAL ANALYSIS OF CARBON DIOXIDE 
RESPIRATION IN MAN* 


ARTHUR B. CHILTON AND RaLpHy W. Stacy 


COMMANDER, Civit ENGINEER Corps, U.S. Navy AND 
DEPARTMENT OF PHysroLocy, Onto STtaTE UNIVERSITY 


By means of physico-mathematical models, formulas are obtained for the 
purpose of making analyses of external carbon dioxide respiration in man under 
conditions of metabolic equilibrium. The methods evolved are flexible, permit- 
ting study of a wide variety of physiological assumptions. Sample calculations re- 
late various factors in the process of external respiration, and show good agree- 
ment with experimental data. A quantitative description of the effects of various 
factors on rate of carbon dioxide output and alveolar tension is given. In particu- 
lar, this theory gives a quantitative prediction of the fluctuation in alveolar car- 
bon dioxide tension over the period of a single respiratory cycle. 


I. Introduction. Until recently physiologists have been unable to ob- 
tain instantaneous values of gaseous tensions within the alveolar spaces 
of the Jungs and associated blood vessels because of sampling and instru- 
mental difficulties. Recent improvements in techniques and instruments 
(Hunter, et al., 1948; Stacy, et al., 1948; Hitchcock and Stacy, 1950) now 
enable better understanding of dynamic changes occurring during this 
cycle. It seems appropriate that a broad theoretica] attack be made on 
the subject to aid in understanding the processes involved. 

This paper initiates such an attack by presenting a logical analysis, 
semi-empirical in nature, of the kinetics of carbon dioxide respiration in 
the lungs under conditions of metabolic equilibrium. To do so, physico- 
mathematical models are set up to represent the respiratory system. 
Actual conditions have been simplified to permit mathematical analysis. 
This simplification has not destroyed applicability of theoretical results 
to practical goals. 

To avoid misunderstanding, we wish to emphasize that this paper does 
not purport to solve the problems of CO: respiration completely. The 
phrase “metabolic equilibrium,’’ which we used above, refers to a condi- 
tion of the human organism in which all the factors involved in respiration 


* The work described in this paper was performed in part under contract between the Office 
of Naval Research and the Ohio State University Research Foundation. 
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have reached a level of adjustment with each other (and with the external 
environment) which permits the breathing cycle to be completely repeti- 
tive in nature. In our analysis, therefore, we may consider the physical 
factors of respiration rate, cardiac output, carbon dioxide level of venous 
blood, et cetera, to be established and constant. More extensive analysis, 
wherein these values are not in an equilibrium condition, requires a treat- 
ment of their inter-relationships and is left as a project for the future. 

Under these conditions, it is proposed herein to use the models to (a) 
analyze carbon dioxide tensions assuming constant alveolar volume, 
(b) generalize to cases of inspiration and expiration, (c) determine cyclic 
changes in tension within the alveoli, and (d) apply formulas obtained to 
study the relative importance of various factors in alveolar carbon dioxide 
tension and output. 

II. Analysis of CO, Flow in Model (Volume Constant). Figure 1 indi- 
cates schematically the physical relationships which are considered to 
exist in the alveoli and associated capillaries. This model accounts for all 
gaseous space in the lungs.-Actually, a simpler model may be used for 
analysis. The space marked Z (since it is neither ventilated nor perfused 
with capillaries) makes no contribution to the problem; space Y is small 
and serves only to cut down slightly the effective blood volume flow. 
Space U is ignored at present for the sake of simplicity but will be included 
later in an extension of the methods derived by the simpler model. Thus, 
the model shown in Figure 2 will be used for the basic analysis. 

In obtaining a suitable basis for analysis, the following additional as- 
sumptions are deemed desirable: 

(A) The compartmentation of lung volume into alveoli may be ignored. 
The CO, tension used represents the average value for all alveoli into 
which there is an active CO, flow. 

(B) The length of blood flow path is considered to be the same for each 
blood volume element. 

(C) Blood is considered a homogeneous fluid with respect to its ability 
to hold CO; both physically and chemically. 

(D) The total cross-sectional area of the capillary path is considered 
constant throughout the active portion of the path. (This assumption is 
convenient, though not strictly necessary.) 

(E) The capillary channel (actually consisting of a very fine network) 
is of such a shape that CO; diffusion lengthwise along the path is slow com- 
pared to CO; diffusion out into the alveolar space and therefore may be 
disregarded. 


(F) The individual alveoli are so small and the walls so interlaced with 
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capillaries that no significant CO, tension gradients exist within the alveo- 
lar space at any instant of time. At least, any difference may adjust itself 
in an insignificant amount of time. 

(G) Variations in the rate of blood flow are small and are rapid enough 
so that a constant mean value may be assumed. 

(H) Gaseous interchange between the alveoli and the “dead space” is 
negligible, except during inspiration and expiration. That is, free gaseous 
diffusion between dead and active space is inconsequential. 


Ficure 1. Schematic diagram of the physical relationships in the lungs. Space U repre- 
sents that portion of the lungs which is ventilated but not perfused with capillaries; V is that 
portion which is both ventilated and perfused; Y is perfused but not ventilated; Z is neither 
perfused nor ventilated. The process of breathing is depicted as a movement of a piston 
(hatched) within a cylinder. ‘Dead space” indicated represents physiological dead space. 


FOR CALCULATIONS 
— USING CONSTANT 
VOLUME, THIS 
PISTON IS ASSUMED 
STATIONARY. 


Ficure 2. Simplified diagram of physical relationships 
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(I) Temperature and pressure variations are negligible. 

Other assumptions will be stated or implied as required in the progress - 
of the analysis. 

Let us now define p as the partial pressure of CO, in the alveolar space 
having volume, V. Let M be the molar quantity of CO, in space V. The 
gas constant is R and T is the absolute temperature. Then from the ideal 
gas law it may be stated that 


pV =MRT. (1) 


We also need a relationship between c, the concentration of CO, in the 
blood, and the corresponding CO, tension, g. We assume one of the fol- 
lowing form, where a and 0 are constants: 


c=aq+. (2) 
Under general conditions we cannot state that ¢ is a function of g only, 
for the CO, dissociation curves, shown in Figure 3, indicate that the 
oxygen tension of the blood is also a variable having some effect on c. How- 
ever, the conditions which exist in the capillaries of the alveoli are reason- 
ably consistent with regard to the oxygen tension, and this fact enabled 
early workers to draw a line from point @ to point @ representing a 
“physiological dissociation’’ curve for CO:. The exact shape of this curve 
is not known, but it is most probably flat as indicated. The assumption 


VOLUME PER CENT CO5 


Ficure 3. A portion of the CO: dissociation curve (modified from Peters and Van Slyke) 
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made herein is that the physiological dissociation curve is close enough to 
a straight line to enable us to choose a relationship in the form of equa- 
tion (2) without introducing an appreciable error in results. This approach 
is not as logical as may be desired, but it provides a convenient way to get 
away from the only alternate—determination of the complex inter-rela- 
tionships between chemical-kinetic and physical] factors involved. 

In solving the problem with constant volume we must define the fol- 
lowing quantities: v is the velocity of blood flow through the capillary 
channel; Z is the length of the channel; and B is the cross-sectional area 
of the channel. Subscripts ; and 2 refer to the conditions of the blood at 
the entrance and exit, respectively, of the capillary channel; é refers to 
the average value of c throughout the capillary channel at any instant 
of time. The subscript » refers to known initial conditions, when time (é) 
is zero. It is apparent that c and q are functions of x and ft; whereas p 
and MM are functions of ¢ alone. 

Physical conditions prescribe that p < q, and therefore CO, is passing 
from the blood to the alveolar space at all times. We see that the rate of 
increase of the amount of CO: in space V is equal to the rate at which CO, 
enters the system at point @ minus the rate at which it leaves the system 
at point @), and also minus the rate at which the CO, content of the al- 
veolar blood is itself being built up. This gives us then: 


(3) 


where 
L 
e=7f cedx. (4) 
0 


To obtain a general analytical expression for c, and thence for €, is quite 
difficult. We may make an approximation, however, which can be ex- 
pected to give a negligible error in most practical cases. It is agreed by 
physiologists that, as a blood volume element passes through the alveo- 
lar capillary path, oxygen concentration comes practically into equilibri- 
um with the alveolar pOs. It is also known that carbon dioxide perfuses 
tissue much faster than oxygen. Hence, we can normally expect that 
CO, tensions in the alveolar blood and the alveolar gas exist almost in 
equilibrium for most of the path length, L. The relationship is indicated 
by Figure 4. (Note: The authors have carried out a more thorough anal- 
ysis of how CO; concentration in the capillaries may be expected to vary 
with respect to x and ¢; but, due to limitation of space, the more qualita- 
tive reasoning above is considered sufficient for the purpose of this paper.) 
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We may say then: 
= 02; (5) 


Gop. (6) 


C2 


Using equations (1), (2), (5), and (6) to evaluate ¢ as defined in equation 


(4), we find that: 


ea M+ 5. (7) 


This expression is substituted into equation (3). It is shown in section III 
that the last term in equation (3), which contains ¢, contributes to dM/di 


x 


Ficure 4. Variation of COz concentration of blood in alveolar capillaries at any given 
instant of time. 


only to a minor extent (on the order of 12 per cent); and this fact also 
minimizes the effect of the error obtained by use of the approximate solu- 
tion for ¢ above. We obtain then: 


dM BVv(a-— d) BakT v 


“di. <VR BLED: VL BIRG ae (8) 
Let: 
Sheets Tek (9) 
= wage Vand 4t 5 
RTa(Bv) ' v 
Then the solution of the differential equation (8) may be found to be: 
—b)V 
M = Myo + SOV y _ gry, (10) 


Again using equations (1) and (2), we may transform the above equa- 
tion to: 


| erat Ere MBAR SS 6 Boe oy 11 (11) 

Ill. Matters in Clarification. Before proceeding from a discussion of re- 
sults under constant volume conditions to those under conditions of vari- 
able volume, it would be well to observe certain facts which bear on both 
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aspects. It may be noted that assumptions and approximations have been 
made on the basis of the numerical size of certain physiological constants 
and variables. Let us, therefore, list some numerical data (found in the 
literature) obtained from experimental and theoretical considerations. In 
attempting to indicate what is “normal’’ we appreciate the fact that a 
wide variation in individual values is likely, due to differences among ex- 
perimental subjects and in experimental techniques. The data included 
herein are selected so as to use a reasonable average, give greater weight 
to more recent findings, and “‘round off” figures slightly for simplicity in 
calculation. 

(A) L/v represents the length of time required for an element of blood 
to traverse the alveolar capillary. The analyses of F. J. W. Roughton 
(1945) indicate that this value under resting conditions is about 0.7 sec- 
onds. For exercise conditions we feel that 0.4 seconds is approximately 
correct. 

(B) Let us consider total lung volume at inspiration as 3440 cc. and 
at expiration as 2890 cc., giving a tidal volume of 550 cc. (at body tem- 
perature and 100 per cent humidity). If we take dead space volume as 
165 cc., then after inspiration V equals 3275 and after expiration V equals 
2725. This gives a median value of V as 3000 cc. In exercise we will as- 
sume V varies from 3500 cc. to 2000 cc., with a dead space of 250 cc. 

(C) From Figure 3 we may arrive at the constants for equation (2). A 
straight line between points @ and (@) is represented by the formula 


Vol.% = 20.8 +2pCO:. Gi29 


In the units useful herein, a is easily computed to be about 32 X 107° 
moles/cc. mm. Hg. 

(D) Under normal, steady-state conditions the pCO; of blood entering 
the alveolar capillaries is a constant, so that we take gq, as equal to 44 
mm. Hg. 

(E) Body temperature used is 310° Kelvin. 

(F) In units used herein, the ideal gas constant F is 62,400 mm. Hg. — 
cc./moles — ° K. 

(G) Assuming that almost all the cardiac output passes through the 
capillaries of the active alveoli, we can see that Bu equals about 5000/60, 
or 83.3 cc./sec. In exercise, we assume 250 cc./sec. 

(H) Roughton’s value for the volume of blood in capillaries at any one 
time (at rest) is around 60 cc. This represents BL. Note that BL equals 
Bv-(L/v), and this relationship checks very well with data taken above. 
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Using the above data we can evaluate equation (9) to find that at rest: 


(sec) . (13) 


We are now in position to introduce and utilize a concept called herein 
“effective volume.’’ It has been noted that in the system represented by 
Figure (2) the CO, is contained not only in the alveolar volume (ignore 
dead space in this regard) but also within the blood of the alveolar capil- 
laries. As has also been noted, an element of blood during most of its 
travel down the capillary channel has a CO; tension practically in equi- 
librium with pCO: of the alveolar air (see Fig. 4). Hence, any change in 
alveolar pCO», as long as it does not occur too rapidly, has concomitantly 
a similar change in alveolar blood tension. 

In certain calculations it will be found that matters will be simplified 
if we replace capillary blood by an equivalent volume of gas, representa- 
tive of the blood in its response to pCO, changes. The total of actual 
alveolar volume plus this hypothetical volume is the effective volume, and 
may be evaluated as follows: 


Moles of CO, in alveolar capillaries = cBL . (14) 


Since 6 moles may be considered mathematically as if permanently bound, 
therefore: 


Available moles of CO: in alv. cap. = (c — 6) -BL (15) 
= @pBL. 
If this CO, were in a gaseous phase, equation (1) would apply and the 
product of p and hypothetical volume would be equal to 
(apBL) -RT . (16) 
Therefore, the hypothetical volume is equal to 
aBL-RT . (17) 
Then, calling the “effective volume” V’, we obtain: 
V’= V + aBLRT . (18) 
For the normal individual at rest the above equation may be evaluated 
by substituting numerical data given previously, in which case we have 
V=V+371 (19) 


with V being around 3000, or about 8 times as large as the hypothetical 
volume. This enables us to justify the statement made following equation 
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(7). If the last term of equation (3) had been ignored, we would in effect 
have substituted the “effective volume”’ for the actual gas volume in our 
calculations—an error of about 12 per cent. 

IV. Analysis of CO: Flow in Model (Variable Volume). Having made an 
analysis under conditions of constant volume, we may now use the results 
obtained in order to generalize to conditions where the alveolar volume, 
V, varies with time. 

(A) Variable Volume without Dilution, as during Expiration. Let us first 
reduce equation (11) to differential form, and by taking cognizance of its 
physical significance we can generalize from the constant-volume case, 
for which equation (11) is valid, to the case of variable volume. If we let 


i A Peay (20) 
and put the equation in differential form, we have 


dy 
Je n0y. (21) 


This indicates that whenever a lack of equilibrium exists between CO, 
tensions of a blood volume element and alveolar air, there is a tendency to 
reduce this difference which is proportional to said difference. The constant 
of proportionality, v, represents the effect of the physical and chemical 
properties of the system and is a function of V as indicated by equation 
(9), derived on the basis of V being constant. The generalization arises 
when we state that the tendency to reduce disequilibrium depends on the 
instantaneous state of the system, in spite of the fact that alveolar volume 
may be changing with time. This permits us to use equation (21) in the 
case of variable volume, provided we express V in the equation as a func- 
tion of time. (Note: This is not strictly true. Analysis by the authors, not 
included herein, has indicated that volume changes tend to produce 
transient effects, which are, however, small and persist only a short time. 
These are therefore ignored.) 

Equation (21) may be solved, then, with the knowledge that » is a func- 
tion of time. The solution is easily seen to be: 


eC nee (22) 
If V is a linear function of time (we shall so consider it in our analysis), 
then equation (9) can be written thus: 


: (23) 


eS s+ wl’ 
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where : 
(2 eats (24) 


ids (25) 


Introducing this relation into equation (22) and solving for y, we obtain: 


(i 


By: = (s+ wit) 1/w* (2 6) 
If we apply the initial condition that y = yo, when ¢ = 0, we arrive at the 
equation: 
S 1/w 
y=» (5) . (24) 
Introducing (27) into (20) and putting vo = gi — fo we obtain: 
RY 1/w 
p= a- (m—p)(— 5) (28) 


(B) Variable Volume with Dilution, as during Inspiration. Inspiration 
is a more complex case because, in addition to the effects previously 
noted, dilution with outside air (so-called tracheal air) or dead space air 
must be taken into account. For this analysis, assume: (1) that pCO, of 
the diluting gas (designated as p’) is constant; and that (2) dilution oc- 
curs slowly enough so that CO; tension throughout the alveolar system is 
essentially in equilibrium. The second assumption means that the capil- 
lary blood, because of rapid diffusion of CO, through the capillary walls, 
feels the effect of dilution as much as does the alveolar air. We therefore 
use the “effective volume’”’ concept. 

If to the volume V’, in which pCO is equal to p, we add a volume AV of 
diluting gas, in which pCO; is equal to p’, then we have increased the re- 
sulting pCO, of the system by an amount equal to Ap. By a simple dilu- 
tion formula we have 


PDE ee yes (29) 
Then, 
Ap _ _b'-? 
AV . V'+AV" ae! 
Taking AV to its zero limit we obtain: 
dp _p'—p 


riattad x. (31) 


CARBON DIOXIDE RESPIRATION IN MAN 11 


Since V’ is a linear function of ¢, it may be expressed as: 


V’=g+ hl, (32) 
where 
g= V,+ aBLRT, (S'S) 
dV 
ks ero @ 4) 
Then we have 
dV=dV'=hedt:; (35) 
therefore, 
= = Paes Pe (36) 
TS 


Adding together the total effects of dilution and CO, passage from the 
blood, expressed by equations (36) and (21), respectively, we obtain the 
following equation: 


A een ey ] 
t Ley, 
h 
(37) 
Se? 
g s+wl 
ae J 


The above equation is made easy to solve because the denominators of the 
fractions on the right-hand side are in a simple ratio relation. That is, we 
can introduce a constant k so that 


E+i=k(stwi). (38) 


In this equation k is defined either as 1/w or as g/hs. By substitution 
of the more fundamental symbols it is easily shown that from either 
definition 


RTaBv 
ripe (39) 
tdi 
We may say then 
ap ish? qi—p (40) 


i k(st+twt) | stwt’ 


The variables are easily separated, and the solution is found to be: 


[Fta-e(itg)] =O (ten. (41) 
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The initial condition to be taken is that p = po when / = 0, so that 


/ —k/k+1 
- ii Poe m 
C= (42) 
sl/w 
Introducing (42) into (41) and noting that k-w = 1, we obtain: 
t —(k+1) 
(otha) —L +h) — hm (2+01-4™ 

p= Fo ——. (43) 

As a special case, if p’ = 0 (as when CO.-free air is inspired), 

wt —(k+1) 
kqi— [kg — po (R+1)] (eee 

p= : (44) 


Raat 


V. Applications of Derived Formulas to Respiratory Problems. The formu- 
las given by equations (9), (11), (28), (43), and (44) provide the means 
of computing the changing alveolar pCO, and the cyclic CO: output of 
the lungs under equilibrium conditions. Let us apply them to the “‘nor- 
mal” data which have been previously listed. Figure 5 shows the respira- 
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Frcure 5. Cyclic variation of alveolar pCO» and alveolar volume in “normal” respiration. 
Some pause after expiration has been assumed in the calculations, so that the figure does not 
indicate the actual agreement between the actual volume curve and the assumed approximate 
volume curve. 
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tory pattern which is consistent with these data. For computational pur- 
poses, the breathing pattern is approximated rather closely by a series of 
straight lines. The computations are made separately for each section of 
the pattern, which can be called phases a-d, etc., to correspond to the 
lettered points on the time axis. Phase c-d represents the first portion of 
inspiration, during which dead space air is being re-inspired. Point d 
therefore represents the instant when fresh air from without is first drawn 
into the alveoli. 


In phase a-d, the expiration phase, 

V = 3275 — 3671, (45) 

and, by use of equations (9) and (28), we can find 
pp =9.02+0.795p, . (46) 
It is necessary to compute the average value of p for the last 0.45 seconds 
of phase a—0. This is the time necessary for expulsion of 165 cc. of air from 
the alveoli, which portion does not escape to the outside atmosphere but 
remains in the dead space to be mixed and re-introduced to the alveolar 
space during the first 0.45 seconds of the next inspiration. This average 
is obtained by integrating equation (28) over the last 0.45 seconds of the 


expiration and dividing by the time increment. This will give the following 
value of p’ for the phase c—d: 


p',= 7.79 40.8239, . (47) 


Phase b-c is the resting phase, at constant volume, and is computed simply 
by use of equations (9) and (11). We obtain: 
b: = 3.5 24-0:920p,)8 (48) 
In phase c-d, the first part of inspiration, 
Vee G7 ie (49) 
and, by use of equations (9), (43), and (47), it is found that 
pa = 3.45 +0.0406f, + 0.8819,. (50) 


In phase d-a, the latter part of inspiration, 
V = 2890+ 3671, (51) 


and, by use of equations (9) and (44), one determines that: 
pa=6.06+0.764p, . (3:2) 


Equations (46), (48), (50), and (52) constitute a set of four equations 


14 ARTHUR B. CHILTON AND RALPH W. STACY 


in four unknowns which can be solved simultaneously to give us: 


ae 


Pe eole3Z | 
eg ckig aap) 


These values are plotted in Figure 5, below the respiration pattern curve. 
The alveolar air which is expelled and lost to the outer atmosphere is 
that portion eliminated from the alveolar space during the first 1.05 sec- 
onds of phase a—b. Let W be the volume of air expelled from the alveolar 

space from beginning of expiration to any arbitrary time. Then, 
W = 367t. (54) 


In the expelled air the amount of CO; in a volume increment, dW, can 
be called dM. 


(55) 


By integrating and substituting proper numerical data, we find that 


M= moles/cycle \ 


vides 
1410 
3279 cofgnin( STP) sl 


(56) 


This figure may be checked by calculating the amount of CO, dumped 
according to the dissociation curve (Fig. 3). From the data calculated we 
judge that the average value of p is about 36.3. If the blood enters with a 
pCO; of 44 mm. and leaves on the average with a pCO; of 36.3 mm., the 
pCO; has decreased 7.7 mm. This represents a difference in volume per 
cent of 5.5. If 5000 cc. of blood pass through the active alveoli per minute, 
the amount of CO, emitted is (5000/100) X 5.5, or 275 cc./min. This 
checks quite closely with the result found by the previous method. 

It may be noted at this point that the compound curve in Figure 5, 
from pa to pa, resembles a simple exponentially rising curve and suggests 
that perhaps adequate results could be obtained by use of a single simple 
formula for that portion rather than several more complex ones developed 
herein. As a matter of fact, equations (9) and (11) will serve very well pro- 
vided that a median value of V is used for the determination of the equa- 
tion of the curve. By employing this method results were obtained as 
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follows: 
Pa = 34.90 (assumed as basis) 


57 
b= 3730 OM 


Comparison with the results, listed as equation (53), indicates the accu- 
racy that may be obtained with this simplified method. Such accuracy of 
course could not be expected if alveolar volume should vary widely over 
the period of the cycle. 

Apparently, the shape of the expiration portion of the respiration vol- 
ume curve (Fig. 5) makes little difference since the rise in pCQ, in the 
alveoli depends primarily upon the median volume in the cycle. 

All the previous results and calculations have been based on what is 
called the “normal” individual’s characteristics while at rest. There are 
wide physiological variations from these characteristics under equilibrium 
conditions of exercise, of varying external environment, or when pathology 
exists. It is of interest therefore to repeat the previous calculations for 
different conditions. Table I indicates the variations in data used and in 


TABLE I 


THEORETICAL ANALYSIS OF THE INFLUENCE OF CHANGES IN SOME PHYSIOLOGICAL 
FACTORS ON CO: OUTPUT AND ALVEOLAR pCO, 


CHANGE IN PHYSIOLOGICAL FACTORS CHANGE IN ComPuTED RESULTS 
eas 5 Median Alv. | Cyclic pCO» 
Altered Factor Amount ene pCOz Variation 
(mm. Hg.) (mm. Hg.) 
1...| Slope of CO2 Diss. Curve +10% + 1% +0.5 —0.1 
2...| Av. Alveolar Volume — 8% 0 0 +0.2 
3...| Venous CO2 Tension +1 mm. Hg. + 2% +0.8 +0.1 
4...| Total blood flow rate +10% + 1% +0.5 —0.1 
5...} Blood Time in Alv. Cap. +14% 0 —0.1 —0.1 
6...| Breathing rate +75% +55% —4.3 —0.5 
Jee Alp lodenate EXercisCces.) alla eeeyecia cor: = OAQU Ai wilt anes +5.7 
8...| Tidal Volume — 9% —11% +0.8 —0.3 
9...| Dead Space — 9% + 33% —0.3 0 
10...| Inact. Part of Alv. Vol. 9% — 6% +0.4 0 


results obtained. In this table, data are given relative to the basic solution 
for “normal” conditions. With one exception, as noted below, the calcu- 
lations follow almost precisely the pattern presented above. 

It will be noticed that for Case 10 of Table I, 9 per cent of the alveolar 
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volume is considered as inactive in the calculations. This means that we 
have used a model more complicated than that of Figure 2, for we have 
taken into account both spaces U and V of Figure 1. The procedure hither- 
to described is complicated by the following facts: (1) during expiration 
part of the expired air comes from the inactive as well as the active vol- 
ume; (2) dead-space air which is re-inhaled is made up of air from U as 
well as from V, in proportion to the relative volumes of these spaces; 
(3) dilution by inhalation affects the air of both spaces, U and V. In spite 
of these complications the calculation is straight-forward and need not 
be explicitly detailed. Results in this case are: 


Poa sa30 

Paes 128 

pe=S1eso 

pa= 38.20 (58) 
fa= fo= fe= 10.29 

fo= 1 .6n 


M, = 257 cc./min. (STP) 


(f refers to pCO» for space U, as p does for space V). 

From the results of the analyses, indicated by Figure 5 and Table I, 
several important facts may be discerned with regard to external CO, 
respiration. 

(A) The result of the basic case—taken as “normality’”’—gives a value 
of CO, output which is reasonably close to that value generally accepted 
as average for human beings. The classic value is about 20 per cent lower, 
but considering the vagueness of some of the basic data, this may be con- 
sidered good agreement. It might be noted that Case 10 of Table I, which 
is probably closer to actual physiological conditions than the basic solu- 
tion, gives a value of CO, output which is lower and closer to the tradi- 
tionally accepted value. 

(B) It is interesting to note the relative importance of various physio- 
logical factors in their ability to affect directly average pCO; of the alveoli, 
cyclic variation in pCO, and total CO, output. Of particular interest is the 
fact that a change in average alveolar volume, or a change in the length of 
time a blood element takes to pass through the alveolar capillaries, has 
little effect on any of the results. Also, we note that an increase in CO, 
tension of the incoming blood does raise somewhat the alveolar pCO, but 
has much less effect on CO, output than one might expect. (This does not 
take into account the effect a change in one factor may have on another. 
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See paragraph below.) On the other hand, an increase in breathing rate 
has a very definite effect on all results. Exercise involves a radical change 
in many physiological conditions and has a great effect on alveolar COp 
respiration. 

A warning must be inserted here. A change in one of the physiological 
factors in a given individual involves more than the direct results indi- 
cated above. If the total picture is desired, it is necessary to take into ac- 
count other physiological adjustments via reflex pathways and chemical 
effects on the respiratory center. These factors have been discussed by 
J.S. Gray (1950). For example, should the cardiac output increase 10 per 
cent, the resulting 0.5 mm. pCO: increase of alveolar air (hence, of ar- 
terial blood) would cause respiratory center stimulation so that CO, out- 
put would increase. Future studies which do take into account such re- 
adjustments among all the respiratory factors are indicated. 

(C) On the assumption that blood as it leaves the alveolar capillaries 
is practically in equilibrium with alveolar pCOz, we see that the blood 
normally undergoes a decrease of about 7.5 mm. Hg. in CO, tension in 
its passage through the lungs. Or, if one prefers to consider the model pro- 
posed by Case 10, this drop becomes about 7.0 mm. Hg. Such a value is 
very close to experimentally determined results. 

(D) One of the most interesting predictions of these calculations is the 
cyclic variation in alveolar pCO:. This fact is one which has recently been 
verified experimentally, though not yet quantitatively measured in man. 
Our calculations show a fluctuation of 2.5 to 3.0 mm. Hg. during the re- 
spiratory cycle. Actually, the sharp breaks in the pCO, curve shown in 
Figure 5 are most likely rounded off. We can estimate that actual experi- 
mental values (should they become available) would be somewhere be- 
tween 2.0 and 2.5 mm. Hg. Note should be made of the large variation 
predicted for exercise conditions. This should indicate to workers an ex- 
perimental condition more amenable to measurement than the normal. 

In addition to providing general Jight on external CO, respiration, as 
indicated above, the results tabulated may be used to solve many concrete 
problems. In general two kinds of problems may be attacked by use of 
Table I. The first type is the analysis of respiratory conditions in an indi- 
vidual whose physiological factors in equilibrium differ from those as- 
sumed by us as “‘normal.’’ The second type of problem is the study of the 
results of factor changes in a single individual. An example of each type of 
problem follows. 

Problem (1). If an individual in a state of metabolic equilibrium exhibits 
all the characteristics listed herein as ‘“normal,’’ except that because of 
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individual variation in the ability of blood to carry CO, the slope of the 
dissociation curve is 20 per cent less, what difference from the “normal” 
might be expected? 

Answer. From the tabulated data, the only difference we can sée is that 
the average alveolar pCO: would be about 1 mm. Hg. lower than that in- 
dicated for Case 1, Table I (normality). This difference assumes that con- 
comitant with the change in slope of the dissociation curve, there is a dif- 
ference of sensitivity of the respiratory center such that alveolar ventila- 
tion is maintained at the level used in the original calculation. 

Problem (2). If an individual, after attaining equilibrium at rest, volun- 
tarily hyperventilates by increasing the depth of respiration 18 per cent 
what is the effect on the results? 

Answer. Such an increase in tidal volume is seen to cause an increased 
CO; output of about 22 per cent, along with a decrease of 1.6 mm. Hg. for 
alveolar pCO, and an increased cyclic variation in same. If the individual 
remained in the original resting condition, physiological changes would oc- 
cur, and his body would not continue the increased CO; output. The nature 
of the compensatory changes is not predicted herein. However, two 
changes might be expected: a reduction in venous CO, tension and a de- 
creased breathing rate. A drop in CO; tension of entering blood of 55 mm. 
and a reduction in breathing rate of 15 per cent would approximately re- 
store the respiratory CO, output to equilibrium with the metabolic output. 

Thanks are due to Professor H. D. Landahl for a suggestion which pro- 
vided a neat and shortened analysis of the basic problem. 
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The conditions are examined under which an approximate relation between 
the radius of the tube and the distance along the axis, as obtained by N. Rashev- 
sky (1945), is consistent with the assumptions made in the solution. These condi- 
tions are reduced to relations between two dimensionless parameters of the sys- 
tem. Second approximations are found for three different ranges of values of these 
parameters. 


Whenever a viscous fluid flows through a tube, there must be a drop in 
fluid pressure along the tube to overcome the frictional resistance at the 
walls. As a result of this pressure gradient and the elasticity of the tube 
wall, the radius of the tube will decrease with distance along the axis. 
Hence the boundary condition of vanishing velocity at the wall has to be 
applied at an unknown boundary. 

This problem was studied by N. Rashevsky (1945) who found a solu- 
tion to a first degree of approximation and then used this solution to ob- 
tain a second approximation. We shall here consider the same problem 
in more detail. 

Fundamental Equations. We use a cylindrical coordinate system 7, 0, 2 in 
which ¢ is the radius and z is the distance along the axis of the tube. The 
cross-section of the tube is assumed to be circular; hence we assume axial 
symmetry, so that the velocity component v, and its derivatives vanish. 
We then have the following fundamental equations. 


Equations of Motion 
i) Equation of continuity: 


1 drv Ov 
a4 r B=) 1 
PAO? Tay ’ (1) 


where 2, and 2, are the radial and axial velocity components respectively. 
We can also write an integral form of the continuity equation, viz., that 
the total mass flow across any tube section of radius R(z) is constant and 
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equal to, say, Q: 


Xe 


R(z) 
2np f rv,dr=Q. (2) 


ii) Navier-Stokes equations for axially-symmetric, steady, laminar flow 
of an incompressible fluid: 


0 0; TT 1 ~3, (G+ ico Us 0? Ve 9) 3 
Urge Te a nag Ney ce (3) 
On, Voy 1 0p Ot, . IL OD: a) 4 
cha yan mer CE ae ara eoe Ce 


where # is the pressure (it is convenient to take p as the difference between 
the fluid pressure and the constant pressure prevailing outside the tube); 
pis the density, a constant; v is the kinematic viscosity = u/p, where pis 
the coefficient of viscosity. 
iii) Force per Unit Area. The normal compressive force on a unit area 
perpendicular to a radius is given by: : 
U; 
Emp (1,2) = p— 2u S—. 


(S) 


Equation of Elasticity 
If R varies very slowly with z and if we neglect effects at the ends of the 
tube, then the forces in the tube will have no component parallel to the 
axis. Then, if we assume a linear stress-strain relation, the following equa- 
tion can be derived, provided 6, the wall thickness, is very small compared 


with R: 
1 1 


Er (R,2) =E3( ge —R), (6) 
where E is Young’s modulus, R, is the unstretched radius of the tube 
(subsequently assumed to be constant), and E,, represents the difference 
between the normal compressive stress on the interior of the tube and the 
external pressure. 

It should be noted that in the derivation of this equation we neglect 
the effect on the radius R(z) of the fluid friction on the interior of the 
tube, and that, due to the slow variation of R(z), we take the component 
of E,, which is normal to the wall to be equal to E,, itself. 

If the strain (R — R,)/R. is so large that the stress-strain relation 
cannot be assumed to be linear, but the variation in R over the length of 
tube considered is small, so that the range of strain along the tube is small 
enough, then (6) can be replaced by 


E,, (R, 2) — 2. =B.8(-Ray), (7) 
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where #, is the constant difference between internal and external pressure 
when the tube is stretched uniformly under hydrostatic pressure; ps has 
a value close to £,,(R, z) in the actual problem; R, is the uniform radius 
of the tube when subjected to #,; E, is the slope of the stress-strain curve 
at the point where the strain is (R, — Ry)/Ru. 

Equation (7) is the result of assuming a linear variation of stress with 
strain over a small range around (R, — R,)/Rx. 

The following treatment uses (6) but is equally valid for (7), provided 
FE, be substituted for E. 

First Approximation to the Solution. If R(z) were a constant then an 
exact solution of equations (1), (3) and (4) would be the flow commonly 
known as Poiseuille flow, in which 

v= 0 and ng Gh (A RY. 
This flow would occur if the tube were rigid; practically, this means that 
the tube would have to be sufficiently “‘strong”’ so that variations in the 
fluid pressure along the tube would have a negligible effect on the tube 
radius. We shall see later more precisely what is meant by a “strong” tube 
in connection with this problem. 

Following Rashevsky, we now assume that in the case of a non-rigid 
tube we can, to a first approximation, neglect all terms containing 2, or 
its derivatives, i.e., that the tube is strong enough so that R(z), though no 
longer a constant, varies so slowly that the resulting radial velocity and 
terms involving it are negligibly small. We shall investigate a posteriori 
under what conditions this assumption may be expected to yield reason- 
able results. We now proceed without further assumptions. 

From (3) we find: dp/dr = 0, i.e., p = p(z) and hence from (5) it fol- 
lows that: E,,(R, z) = p(z). Then (6) gives us 


dp _EaaR i 
HEARED SL 

and (4) becomes 
1dp a4, Lia j 
pas” eg 2 ae : oe 


Integration of (9) and application of the boundary condition v, = 0 at 
x = R yields: 


Se UA ee ie (10) 


aie 


For any fixed z, (10) gives a parabolic radial distribution for v,. Hence 
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we see that if we neglect terms containing v,, then the velocity profile at 
each section must necessarily be the same as in Poiseuille flow. 
Substitution of (8) into (10) gives us: 
Eé6 dR r2 
aaa peleg ( | (ABE 11 
a 4u dz : 4 Ose 

Finally we apply (2), the continuity condition, and obtain: 

dR 8uQ (12) 


dz Eon pR? 
and 
RiBR peek 
3 SaaS ee, 
where Rj is the radius at z = 0. 

Equation (13) is identical with Rashevsky’s result. The method of its 
derivation, however, differs slightly from Rashevsky’s in that the parabolic 
variation of v, with radius [equation (10)] was shown to be a consequence 
of the assumption concerning the smallness of v, and its derivatives, rather 
than an additional assumption. 

Validity of the Approximation. We can now examine under what cir- 
cumstances the above solution will be consistent with the assumptions 
made provisionally in the section above. 

The exact equation (4) may be replaced by the approximate equation 
(9) provided 


ree eee 
are very small compared with 
2 
: oH 52): 
From (8), (9), and (12) we have 
O70, 1:09, 8 uO 
"Nar vires ow pR* 


Also, using our results, the following orders of magnitude are easily 


verified: 
02 u, 303 
aye eS €: pam) 


Ome _ rey 
& Or =0 rae) 
4 OU, - uQ3 ) 
a8 9 asp 
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Hence we may neglect the non-linear terms, provided 


QV? 
Es pRi<'° (14) 


Or, writing Q = pV R? where V is an average velocity over a cross section, 
we have 


See (15) 


elie (16) 


or 
uV\" (17) 


If we write pVR/u = §, a Reynolds number appropriate to our prob- 
lem (V and R may be taken as the values at z = 0) and pV /E6 = &, 
then conditions (15) and (17) become respectively 

RN<KI (18) 
and 
N2< 1. (19) 

We note that Jt is a dimensionless ratio of a viscous to an elastic force. 
We may regard # and Jt as the two dimensionless parameters of the sys- 
tem. For any fixed 9% the flow will approach Poiseuille flow as It — 0. 
Our provisional assumptions also involved equations (3) and (5). We can 
write OE,,/dz = 0p/dz provided |u(d%,/drdz)| « |dp/dz| and this, 
using our results and equation (1), is found to be true if 2? < 1. 


Finally, using (3), we find 
10p 2%w0dR/(8 °) : 
a ere eA + higher order terms. 


By integrating from 7 to R and using (6) we obtain 


_ns(1_ 1), 8x0 aR or) 


whence 
az R? pra Tp — zt Tp a2? Pe ae 
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The last two terms are negligible, provided y?Q?/E°6’p?R* <1 or 
M2 < 1. Thus (18) and (19) are the conditions which must be satisfied. 

Let us first consider the case of 3 >> 1. If (18) is satisfied then (19) 
will be automatically satisfied, since it is a weaker condition in the pres- 
ent case. Hence, for 9 >> 1, the results are consistent with our assump- 
tions if Jt is so small that RM <1. 

If 9% = O(1) then (18) is true if Jt <1. Hence in this case this is the 
only condition. 

If R <1 and R = O(M) the results and assumptions are consistent 
1B 

We now see that a “‘strong”’ or “practically rigid” tube in connection 
with this problem means that #Yt and MN? are very much less than one. 

Second Approximation. Cases R > 1, RN K land RK = O(1),N «K 1. 
We can obtain better results than the foregoing by taking account of some 
of the terms previously neglected by a method of successive approxima- 
tions: we shall include terms of second order by computing their approx- 
imate value from the first approximation. Calling the terms retained in 
the first approximation terms of order one, the next higher order terms 
in the present case will be those of order RY, i-e., the quadratic terms in 
(4). All other terms are of order 9? which is negligible compared with RY. 
Hence we still have 0p/dr = O and (8) holds. It is important to note that 
the quadratic terms 2,(0v./dr) and v,(0v./0z) are both of order RY and 
that both must therefore be included in the present approximation. 

Using the results of the second section and equation (1) these terms 
can be written as follows: 


Vz _ 8Q? dR 1 Lire 
Sf UP soe Ce ’ Z) - wp? dz ARS R? 5): (223 
Hence (4) in our second approximation becomes: 
1dp oa 10, 
I(r2) = 2 ty ($3 +55"), (23) 


Our first solution resulted in a parabolic radial distribution of v,. This 
was essentially due to the fact that in the approximate equation (9) the 
viscous term was not a function of 7. In (23) this is no longer the case 
and hence we expect a modified velocity profile. Integrating (23), using 
(8), and applying the boundary condition v, = 0 at r = R, we obtain 

E6 dR 20? dR re ra aa 
etre ers) ae da\~ oR * IRI ne ae 


Comparing equation (24) with (11) we see that the second term in (24) 
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is the additional term introduced by taking account of the inertia terms. 
Finally, applying (2), the following equation is obtained for dR/dz: 


Eirp ps, Q 1) dR_ 


8uQ QrpR/ dz >) 


Since dR/dz is negative we see that the effect of the correction term 
Q/2ruR [compared with the corresponding term of (12)] is to cause R to 
decrease more rapidly. Equation (24) shows that for a given R(sz) the cor- 
rection term tends to reduce the speed on the axis and to maintain the 
given volume flow by a sharper gradient near the wall. 

We note that (24) and (25) hold subject to the previously discussed 
restrictions; moreover, since the conditions (14) and (16) involve the vari- 
able R, the equations only hold over a range of z in which R does not be- 
come too small. 

The above case is the one treated by Rashevsky, who obtains an equa- 
tion identical with our (25) except for the fact that the numerical coeffi- 
cient $ of his correction term is replaced in (25) by 4. This discrepancy is 
due to the fact that Rashevsky only takes account of the second of the 
non-linear terms in (4) and, moreover, that he applies some average value 
of this term over the cross-section. This latter step results in the viscous 
term of (23) being independent of as it is in the first approximation and 
hence in the retention of the parabolic velocity profile. 

Second A p proximation. Case  « 1. Here the first approximation would 
hold even if 3 = O(9). The terms to be taken into account in the second 
approximation will depend on the order of N. If Jt « #, then a second 
approximation will hardly be useful since the correction terms will be very 
small. If Jt = O() then, in the second approximation, we must take ac- 
count of terms of order RY and N?, i.e., of equation (3) as well as all terms 
previously neglected in (4) and (5). 

This can be done by a method analogous to that used in the preceding 
section by replacing all unknown terms of second order by their approxi- 
mate values as obtained from the first approximation. 

Using this process, equation (4) becomes 


1(r,8) = 2 Shy (Fa 4e Sten), (26) 


where G(r, z) is the known first approximation of 0°2,/0z’. Also we have 


dp EédR 
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where H(r, z) is a known function which corrects for the fact that 0p/dr 
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is no longer negligible. The function H is determined from equations (3), 
(5), and (6). 

If we introduce expression (27) for 0p/dz into (26), integrate with re- 
spect to y and apply the boundary condition on r = R, then application 
of (2) will yield the new equation for dR/dz. It is: 


Eom p a, oe HO els Caen eke 
ma - R2 =], 2 
SiO 1) Sekona Ri Wainy Iemma 22) 


Because of the very small values of 3t which are likely to be encountered 
in practice, it is to be expected that for most applications Jt < 9; in this 
case the second term in parentheses in (28) becomes negligible compared 
with the third term and the equation reduces to (25). 
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An hypothesis regarding the mechanism of spatial representation in the neural 
centers is formulated in order to explain the audiogyral illusion. Using this hy- 
pothesis and experimental data (Clark and Graybiel, 1949) the time constant of 
the semicircular canals is calculated. The value obtained in this manner agrees 
with that previously calculated (Mayne, 1950) using results from experiments on 
the nystagmic latency of pigeons. 


The Audiogyral Illusion. Clark and Graybiel (loc. cit.) determined the 
effect of angular rotation on the localization of sounds using human sub- 
jects. A mistaken judgment regarding the direction of sound was noted 
during the period immediately following rotation. The phenomenon was 
called an audiogyral illusion. 

The paper by Clark and Graybiel gives a complete description of the 
experimental procedure. In brief, however, the subjects were first rotated 
at a constant velocity in a Barany chair, then stopped suddenly. They 
were then asked to judge the apparent direction of a sound source located 
at a fixed position with respect to themselves and to the chair. The appar- 
ent displacement of the source from its actual position was noted as a 
function of the elapsed time after rotation was stopped. A total of six 
series of tests was made, for clockwise and counterclockwise rotation and 
for three initial positions of the sound with respect to the observer—ten 
degrees to the right, ten degrees to the left, and center. 

A number of control runs were also made to evaluate the ability of the 
subjects to estimate correctly the direction of the sound under static con- 
ditions, i.e., without previous rotation. It was found that the mean error 
for the three positions was only two degrees. 

The results of the test are best summarized in the words of the origi- 
nal paper: 

Clearly defined changes in sound localization occurred following the angular 


decelerations used in this experiment. The changes were consistent for the three 
subjects and showed that after abrupt decelerations, they tended to localize the 
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sound as if it were displaced in the direction of the preceding rotation. Thus, after 
rotation to the left, a sound coming from a source directly in front of the S’s 
would be reported as coming from the left of center. On the other hand the illu- 
sory feelings of rotation would be opposite in direction, i.e., the S would report 
that he felt as if he were rotating to the right. 


An attempt will now be made to interpret the results and to fit these re- 
sults into a theoretical framework of the behavior of the semicircular 
canals and of the neural representation of space. 

The Post Rotational Effects. Since the audiogyral illusion takes place 
during the period immediately following the cessation of a steady rotation, 


Ficure 1. Schematic representation of semicircular canal (Jour. Comp. and Physiol. 
Psychol., 43, No. 4, August, 1950). 


we should consider briefly the behavior of the semicircular canals during 
this particular period. 

In a previous paper (Mayne, 1950) the thesis that the semicircular 
canals issue signals which are proportional to the angular velocity of the 
head was developed. It was shown that this proportionality obtains only 
within certain ranges of amplitude and frequencies corresponding to nor- 
mal body movements. During the post rotational period the semicircular 
canals issue definite velocity signals although the head is stationary. A 
quantitative prediction of these post rotational signals can be made with 
a brief consideration of the mechanism of the semicircular canals. 

Referring to Figure 1, illustrating diagrammatically the structure of the 
semicircular canals, the following equation can be seen to describe the mo- 
tion of the endolymph: 


d*6, d?6 dé 
(Gr Ge) Gt Ke. (1) 


This may be written as: 


dO _@0 , cd0, Ké 
di Gee ee 2) 
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m = mass of the system 


¢ = viscosity of the system 


K = elasticity of the system 
9, = displacement of the walls of the canals in space 
6 


= displacement of the fluid with respect to the walls. 


From the previous paper (Mayne, 1950) we have: 


= 0. 
m 
Lae 
m 
Equation (2) then becomes: 
765 a6 dé 
dp Sayer 200 sa ls 246. (3) 


Figure 2 illustrates the displacement of the fluid, as calculated from 
equation (3), when the body is accelerated abruptly to a constant velocity; 
then, following a period of half a minute or more, is suddenly stopped. 


————-— SENSED VELOCIT 
ACTUAL VELOCIT 


ANGULAR VELOCITY 


Frcure 2. Sensed angular velocity for step velocity input 


This assumed velocity input is represented by the solid lines in the dia- 
gram. The dotted lines represent the displacement of the endolymph and 
the velocity signal corresponding to this displacement. The effect of the 
first term of the second member of the equation can be neglected at the 


scale of the diagram. 
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Calculations indicate that this term has the effect of delaying the abrupt 
rise of the fluid displacement by .015 sec. 

The return of the endolymph and of the cupula to normal position fol- 
lowing abrupt displacement, neglecting as mentioned above the accelera- 
tion term of the equation, is given by 

6= He’, (4) 
where 

§ = displacement of the endolymph from neutral position at a time ¢ 

following abrupt acceleration or deceleration, 


6) = initial displacement of the endolymph at the end of acceleration 


period, 
ee ¥ = the time constant of the semicircular canals, 
¢ = time from the end of abrupt acceleration or deceleration, 


e = base of natural logarithm. 


The curves show that if the body is accelerated abruptly to a constant 
velocity, the canals will give, immediately following the acceleration, a 
signal normally corresponding to this velocity. The signal then decays to 
zero according to equation (4) even though the body is maintained at a 
constant velocity. If rotation is now suddenly stopped the same phenome- 
non takes place, except that the endolymph is displaced in the opposite 
direction. During the return of the fluid and of the cupula to normal posi- 
tion the canals issue velocity signals even though the head is stationary. 

These false signals give rise to conflicting information and the well- 
known post rotational phenomena. They explain the sensation reported 
in the last sentence of the quotation from Clark and Graybiel. They do 
not explain, however, the illusory displacement of the sound source. A 
hypothesis of the cause of displacement must be formulated. 

The Role of the Vestibule in Preset Movements. It is known that certain 
movements of the body can be carried out without the continuous assist- 
ance of external senses. This type of control has been discussed, and it is 
shown (Mayne, 1951) that such movements are carried out in closed loop, 
with preset input in the neural centers. Such movements are said to be of 
the preset or autopilot type, to distinguish them from those movements 
which may be carried out under continuous control of external senses. In 
the case of a preset movement of the hand, the input in the neural centers 
may be matched with signals of the proprioceptive senses for closed loop 
operation. In gross body movements the neural input must be matched 


AUDIOGYRAL ILLUSION 31 


with signals from an absolute spatial reference. It is believed that the 
vestibule supplies such a reference by a process of double integration 
(Mayne, 1950, 1951). The first integration is conceived as taking place in 
the vestibule proper and the second in the neural centers. 

The stability of the movement demands that it be controlled as a func- 
tion of both the error and the velocity of the movement. The error, as in all 
servo-systems, is the difference between actual and desired position. It 
is possible to conceive of a system of spatial representation and of a re- 
lated operation of the vestibule, so that the above requirements for posi- 
tional and velocity signals will be consistent with the known behavior of 
the semicircular canals, with post rotational effects, and the presently 
discussed audiogyral illusion. 

Spatial Representation. Space, and the arrangement of objects about us, 
is probably represented in the neural centers by means of a digital bank 
of on-off elements. For the sake of simplicity let us consider, rather, an 
analog type of representation. Let us suppose that the various points of 
reference about us are merely spotted on a plane, and let us limit our con- 
siderations to a two-dimensional representation. Let us draw on our 
plane, first, two perpendicular axes, representing the front and back, and 
right and left axes of the body. The eyes may be assumed to be fixed with 
respect to this system of axes, and the problem may be further simplified 
by assuming that the head is rigid with respect to the body. 

Now, the necessary information is supplied through the complex process 
of vision to spot various points of reference on the plane. The representa- 
tion is easiest if the spatial arrangement has symmetry about two axes, 
and if these axes coincide with the body axes as represented on the plane. 
This may be why arrangements constructed by men are usually sym- 
metrical about two axes, and why observers naturally orient themselves 
so that their body axes coincide with those of the arrangement. 

The spatial representation may include any number of reference points 
or lines. At any instant we can close our eyes and have sufficient stored 
information to perform any number of tasks: touch the corner of our desk, 
reach for the telephone, stand up and walk to the door of our office with- 
out bumping into a chair, etc. For the purpose of this discussion we can 
reduce the information spotted on the reference plane to a set of two per- 
pendicular axes, one set corresponding to the points of the compass, or to 
the main axes of symmetry of a room where we may be located, and an- 
other set corresponding to the body axes. Our reference plane will then 
have two sets of axes: an absolute set referred to space and a relative 


set referred to the body. 
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If at any time we turn our body with respect to our surroundings, the 
relative axes must be rotated with it, but the absolute system must re- 
main fixed in space. If, for instance, we sit on our revolving desk chair 
with our eyes closed and someone turns us 90 degrees, we still can touch 
our nose and therefore must be aware that it has changed its location in 
space, or that it is in the same position with respect to our body. At the 
same time we can still reach for the telephone on our desk, so that we 
must be aware that it has remained fixed in space, or has changed its loca- 
tion with respect to our body. 

Proper orientation requires, therefore, that the absolute axes be rotated 
with respect to the plane of representation and to the relative axes during 
rotation of the body. This rotation is believed to be effected as a function 
of signals from the semicircular canals. It is evident that there are no 
planes in the neural centers where space may be represented as sketched 
out above. However, an equivalent digital system, where the coordinates 
of a point would be represented by association with two numbers, can be 
easily conceived. 

A desired rotation of the body can then be preset by calling for a certain 
angle between relative and absolute axes as an input. The movement may 
then be controlled by the difference between this input and the actual 
angles between the two sets of axes. This difference would then corre- 
spond to the positional control function. A value proportional to rate 
should be added to this positional function for stable and rapid control. 
Such a rate could be obtained by differentiating the positional informa- 
tion. It appears to be obtained, however, by a more direct method. The 
relative system of axes seems to be shifted ahead by an amount propor- 
tional to the velocity indicated by the semicircular canals, so that the 
error signal would then contain both position and rate terms. 

Computation of the Time Constant of the Semicircular Canals. If the 
foregoing hypothesis is correct, the relative axes of reference of a subject 
are shifted in direct proportion to the velocity indicated by the semicircu- 
lar canals and in the direction of this velocity. A sound referred to this 
system of axes would then appear to be shifted in a direction opposite to 
that of rotation. The shift is not a function of actual velocity, but of the 
signal issued by the semicircular canals, as discussed earlier. In the case 
of the post rotational period, it has been shown that the signal is given by 
the exponential expression (4), and it therefore should be expected that 
the sound shift would follow a similar expression. Further, a plot of the 
shift vs. time should make it possible to compute the time constant of the 
semicircular canals. 
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For the purpose of analysis, the data from the six runs in the study by 
Clark and Graybiel were averaged. There are differences between indi- 
vidual runs which may be significant, and the tests should be repeated 
with more elaborate equipment. But, since each run represents averages 
of a large number of individual tests, there appears to be little possibility 
of interpreting these differences. A closer attention to individual test vari- 
ance in future experimentation may provide further clues to some of the 
details of the system of spatial reference. 


CURVE CALCULATED FROM 
@= 16.9% ¢ - 2.094* 


© EXPERIMENTAL POINTS 


REFERENCE | 


APPARENT DISPLACEMENT OF SOUND SOURCE FROM 
ACTUAL POSITION (DEGREES) 


sd) 5 10 15 20 25 30 35 
TIME FOLLOWING CESSATION OF ROTATION (SECONDS) 


Ficure 3. Experimental and calculated shift of sound during post-rotational period 


As the case may be, the curve of the average disorientation with respect 
to time is approximated very closely by an exponential where 1/7 = 
0.094 1/sec. 

Figure 3 shows the agreement between the experimental and the fitted 
curve. The value of 1/7 = 0.094 computed here for human canals agrees 
reasonably well with that of 0.120 computed in the previous paper, using 
results (Mowrer, 1935) from experiments on the nystagmic latency of 
pigeons. , 

It may be pointed out that objects observed visually do not experience 
the same illusionary shift of position that occurs in the case of sounds. 
This may be explained by the nystagmus of the eyes, which compensate 
for the shift of axes, and in fact overcompensate for such a shift, as shown 
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by other experimental data. This overcompensation is probably for the 
purpose of providing for a reaction-time delay in spotting reference points 
in the neural centers while the head is in motion. It is hoped to present a 
discussion of nystagmus in a later paper. 

Conclusion. The above results lend weight to the following views re- 

garding the semicircular canals and the system of spatial reference. 

1. The semicircular canals are a structure defined by equation (2). 

2. The values of the time constant of the semicircular canals for man 
and pigeons are similar. 

3. The time constant of the semicircular canal has an approximate 
value between 0.094 and 0.120 1/sec. 

4. The semicircular canals give a signal proportional to the angular 
velocity of the head within certain limits of velocities and fre- 
quencies. 

5. The semicircular canals give rate information to stabilize certain 
body movements. 

6. In preset movements rate information combined with positional 
information is obtained by comparison of an absolute system with a 
relative system of axes in the neural centers. The absolute system of 
axes is maintained fixed in space while the relative system is shifted 
ahead in the direction of the rotation, as a function of information 
supplied by the semicircular canals. 

7. This mechanism explains the audiogyral illusion. 
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The spread of excitation in a ‘‘random net” is investigated. It is shown that if 
the thresholds of individual neurons in the net are equal to unity, a positive 
steady state of excitation will be reached equal to y, which previously had been 
computed as the weak connectivity of the net. If, however, the individual 
thresholds are greater than unity, either no positive steady state exists, or two 
such states depending on the magnitude of the axone density. In the latter case 
the smaller of the two steady states is unstable and hence resembles an “‘ignition 
point” of the net. If the initial stimulation (assumed instantaneous) exceeds the 
‘{gnition point,” the excitation of the net eventually assumes the greater steady 
state. 

Possible connections between this model and the phenomenon of the ‘‘pre- 
set” response are discussed. 


Consider a random net Jt such as was described by R. Solomonoff and 
A. Rapoport (1951), i.e., an aggregate of neurons, each emitting a axones 
which synapse randomly on other neurons in the aggregate. Let the thresh- 
old of each neuron be 1. Furthermore, let the postulate of “quantized 
time” hold, that is, firings occur only at multiples of a unit of time, which 
is taken to be the synaptic delay, assumed constant throughout the net. 

Now let a subset of the neurons of Jt be stimulated at ¢ = 0. We wish 
to follow the ‘‘expected”’ course of excitation in Jt. 

If x(é) is the number of neurons firing at the instant ¢ (an integer), then, 
following arguments similar to those developed in the earlier paper (Joc. 
cit.), we obtain the equations 


C0 =w{i-(1-3) |: (1) 
aU D1 exp a (2) 


If a steady state exists, we will denote the number of neurons firing 
steadily by Vy = «(@), and obtain an equation in y identical to equation 
(20) of the paper mentioned above, namely, 

Oct Va SLO (3) 
35 
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This result can also be obtained by writing equation (1) as a differential 
equation, which holds approximately for the case of quantized time: 


dx a pa ae + 
ow (1-exp SE )-2 0, (4) 


Here the right side represents the difference between the expected number 
of neurons firing at ¢ + 1 and the number firing at ¢. But that is the rate 
of change of « per unit time. Hence it is equal to the left side. Separat- 
ing the variables, we obtain 


Ghat. z =p 
N (1—exp aes 


whose solution for ¢ in terms of x is obtained by a quadrature: 


t, (S) 


t= 2 (6) 


The steady state is obtained by setting (4) equal to zero, which again 
leads to 
y= ile". 


where y = «/N at the steady state. If this steady state is actually reached, 
then the net 9 reaches that level of activity regardless of initial conditions. 
The steady state is thus determined by the structure of the net and not 
by the initial stimulation imposed upon it. 

Remark. This model can be somewhat modified by assuming a refrac- 
tory period 6 > 1. Suppose, as an example, 1 < 6 < 2, i.e., no neuron can 
fire at two consecutive units of time, but a neuron may fire every two units 
of time. In that case, equation (5) becomes 


d: — ax 
Giz (N—2)(1-ep S)-«. (7) 


Accordingly, the spread of excitation is given by 


= d 
= (9) 5 s Sx ’ (8) 
‘ (N — 2) (1—exp )-2 


cay 
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and the steady state by the roots of 


7 , — ee ——y 
(N (1 exp V ) CaO (9) 
or 
ye i= : (10) 
Solving equation (3) for a in terms of y, we obtain 
1 1 
while equation (10) gives 
ae iL == Gy 


Equation (12) has meaning only for y < 4, which is intuitively evident, 
since under the conditions of the refractory period above we cannot have 
more than one half of the neurons firing twice in succession. However, it 
can be seen that as a grows to a moderately large value (5 or 6), y ap- 
proaches its asymptotic value of $ very closely. The situation is essentially 
similar to the one without the refractory period, except for the reduced 
limiting value of y as a function of a. 

The “Ignition” Threshold. In the case just considered where the thresh- 
old of each neuron is unity there is no lower limit on the amount of initia] 
stimulation which eventually leads to the steady state activity. Theo- 
retically, even the firing of a single neuron eventually propagates so as to 
bring the net into an active steady state, provided a > 1. This is not sur- 
prising in view of the fact that for a > 1, the number of neurons firing at 
successive times propagates approximately exponentially in the beginning 
of the process. The situation is essentially different if the thresholds are 
greater than unity. 

In the previous case we had for the probability that a neuron receives 
a stimulus at time ¢ + 1 the following expression: 


1 ax(t) — ax 
= Sa mwil= ee 13 
1 (1 x) 1 — exp W (13) 


This result is obtained by the following reasoning. At time ¢, ax stimuli are 
being sent out by the firing neurons. The probability that a given neuron 
is missed by all of them is 

(1-5 eee (14) 


N 


Therefore the probability that it receives at least one of them is given by 
equation (13). However, we can also reason in another way. The proba- 
bilities of receiving no stimulus, exactly one stimulus, exactly two stimuli, 
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etc. are given successively by the terms of the Poisson distribution 


SS (OHAY 
p (0) = exp > — 


— 0X 


ax 
pl) = itv? Ww 


(15) 


akxk — ax 


PAEY ain oP 


The sum of all these probabilities is, of course, unity. Therefore the suc- 
cessive probabilities of receiving at least 1, 2, etc. stimuli are 


P(1) =1-ep = 


P (2) =1—exp—7* 145% 
: (16) 


Pikjo=] — exp Es (+), 
where 
k gi 
By (2), 2 pa 


7=0 


as previously defined in an earlier paper (Rapoport, 1950). We are thus 
led to the generalizaton of expression (13) for the case where the thresh- 
old is k > 1, namely, the probability of receiving at least # stimuli at 
time ¢ + 1, when x(¢) neurons are firing at time #: 


P(h) =1-exp—** (+). (17) 


This leads immediately to a generalization of equation (6) for this case, 
namely, 


dx —ax ax 
GaN [1- exp SB (S)]-«, (18) 


whose formal solution is given by 
ea yh dé 


= 2(0) yy | 1 — exp =H ma (SP) |-§ 


(19) 
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The steady state is given by the roots of 
b= 6 VE (ay) —y=0. (20) 
If 4 = 1, equation (20) solved for y has exactly one positive root for 
1 <a < o,as has been shown (Solomonoff and Rapoport, 1951). It will 
now be shown that if # > 1, the corresponding equation (20) has exactly 


two positive roots, provided a exceeds a certain value, which depends on hf. 
We write equation (20) as 


1 —y = e “EF, (ay) . (21) 
Holding a and / fixed, plot the left and right sides of (21) against y (Fig. 
1). Clearly the abscissae of the intersections of the two curves will de- 
-aZ, 
e E,., (2) 
OR |-% 


ov 


Ficure 1. The line 1 — y and the curves e~*7E;_(ay) for two values of a plotted against y. 
For a sufficiently small value of @ (upper curve) no intersection occurs. For sufficiently large 
values, there are exactly two intersections. 


termine the values of y which will satisfy (21) for a given pair of values 
of a and kh. We note that if # > 1, then 


< e-E,-1 (ay) = —ae—ME,-1 (ay) tae“E,-2 (ay). (2 2) 
At y = 0, 
Ma e-“E,-1 (ay) = —ata=0. (23) 
dy 
Moreover, 
ad? ss 
dy 6 UE, = 1 (ay) = o?¢ “H,_1 (ay) — 2 @e—-VE,—2 (a7) (24) 


+ ae—ME,—3 (ay) .* 


* Consistent with the definition of ;, the value 0 is assigned to Z- (cf. Rapoport, 1950). 
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If we set the right side of (24) equal to zero, we find that since 
Ge 70: 
Fy-1 (ay) — 2Ey-2 (ay) +£i-s (ay) =0, (25) 


which may be rewritten as 


Ga aaa a: - 
(Fn-1 = [i =>) =" (Ej,—2 — Fs) = (h =e 1) fees (h sok 2) ! aE 0 . (26) 


The roots of (26) are at y = 0 and y = (# — 1)/a. This means that 
the graph of e-*’E,:(ay) has at most two inflection points, namely, at 
y = 0 and at y = (h — 1)/a. Consider now the slope e~E,a(ay) at 
vy = (h — 1)/a. Substituting this value into equation (22), we find: 


d 
ar CH= (Ay) \y=Ge/ 
vey, 


= — ae-*DE,_, (h—1) + ae-@ DE,» (hk — 1) (27) 
eee Ch ia riubes. 
Gin V/2r(h—1)’ 


where i > 1. 

Hence a/‘V 27(k— 1) is approximately the greatest absolute value of the 
slope of e~*7E,_1 (ay). If this value does not exceed unity, that curve will 
never cross the line 1 — y, and equation (20) will have no positive root 
(cf. Fig. 1). Therefore we have deduced 


Lemma 1. If kh > 1, no positive steady state exists, unless approxi- 
mately 


a> V2 (h—1). (28) 


Remark. Inequality (28) is a necessary, but not a sufficient, condition on 
the axone density of the net for the existence of a positive steady state. 

Lemma 2.1fh > 1, there exist either no positive steady states or exactly 
two such states.* 

Proof. We note that the right side of equation (21) approaches zero 
monotonically from above. Therefore if that curve has crossed the line 
1 — y at some point where y < 1, it must cross it again at some other 
point where y < 1. The two intersections determine the values of y, 
where dx/dt = 0, i.e., steady states. 

Lemma 3. Let the two positive steady states (if they exist) occur at 


yiand 2 (v1 < yz). Then x/N = y; gives an unstable steady state while 
x/N = 2 gives a stable one. 


* We are disregarding the very special case where tangency occurs. 
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This follows immediately from the inequalities 


dx a6 dx 

— — I) 5 w 

a7 for ips ta apo for NaS: 
ee Osa! a . 
al ; or Wr v2: 


We thus have the following situation. Suppose a and / are fixed in such 
a way that two positive steady states exist. Then if the initial stimulation 
of J is such that «(0) < Ny, the activity will die out. If, however, the 
initial activity exceeds Vy, the activity will increase until steady state 
N72 is reached. If the initial activity exceeds Ny2, it will decrease to that 


STEADY STATE OUTPUT 


INSTANTANEOUS INPUT 


Ficure 2. The all-or-none character of the steady state of excitation in a random net 


steady state. We have here essentially an ‘‘ignition”’ phenomenon, where 
the threshold of ignition acts as the overall threshold of the net. If the 
initial instantaneous stimulation of the net involves a number of neurons 
in excess of that threshold (i.e., Vy;), the net will “ignite” and will remain 
active to the extent of the involvement of Ny2 neurons in that activity. 
The input-output curve of such a net is shown in Figure 2. The input, 
however, must be here understood as an initial instantaneous input and 
not a continuous one as was the case in our earlier papers dealing with 
input-output problems. The output, on the other hand, is the steady state 
output of the net. In terms of these definitions, the net exhibits an ali-or- 
none behavior. 

Relations between the Parameters of the Net. Our treatment has given rise 
to the following parameters, characterizing the net {t: the axone density a; 
the individual neuron threshold /, which is an integer designating the 
number of stimuli which must be ‘‘simultaneously” received by a neuron 
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to elicit a firing; the net threshold 71, designating the fraction of the neu- 
rons which must be initially stimulated for the net to become permanently 
active; and ‘2, the fraction of neurons active in the stable positive steady 
state. We seek relations among these parameters. 

It is easy to verify that if h < h’, then 


e-%E,-1 (ay) < e-%Ey’-1 (¢) (30) 


for all values of y > 0. Therefore, for fixed a, the curve e~”Ey(ay) 
will lie wholly above the curve e~*’£,1(ay), as shown in Figure 3. If 


-ao 
e E,_,@?) 


x 


deat yy 


Frcure 3. The curves ¢-?7E;,_,(ay) for two values of #. The upper curve corresponds to the 
higher value of h. 


y, and ‘yz are the threshold and steady state parameters associated with 
h’, then we shall have 
Yio ay yom (31) 
1 2 


In other words, if the thresholds of the neurons comprising the net Mt are 
uniformly raised, the overall threshold of the net will be raised and the 
steady state activity lowered, as, of course, should be the case. Vice versa, 
a lowering of the individual thresholds will lower the overall net threshold 
and raise the steady state activity. 

We have here a model of what has been called a “preset’’ response 
(Mayne, 1951). Suppose 9 is innervated by other centers, C1, Co,... C; 
in such a way that each neuron of Jt receives an axone from each of the 
C’s. For simplicity of discussion, we still assume that all firings are syn- 
chronized and occur once per unit time. Then if the centers C;.... C; 
are totally activated, and k < h, the stimulation of 0 by the C’s is sublim- 
inal. However, the effect of the continued activation of the C’s is to re- 
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duce the individual thresholds of the neurons of Jt by & units with re- 
spect to stimulation from other sources. Thus the amount of response 
which will occur in 9t as a result of outside stimulation has been “preset” 
by subliminal stimulation from the centers C,.... Cz, and this response 
is again independent (after the steady state has been established) of the 
magnitude of the outside (superthreshold) stimulus. 

We may now generalize this model by supposing that the individual 
thresholds in %t are not uniformly distributed. Or we can suppose that 
the centers C innervate 3t only partially, each center innervating its 
particular set of neurons in Yt. Then a distribution of activity in the cen- 
ters C amounts to a new distribution of individual thresholds in 2. 

Let the fraction f; of the neurons in St have threshold 1, the fraction fe 
have threshold 2, etc. If some values of the threshold are not repre- 
sented in the net, the corresponding /’s can be taken to be zero. Hence 
formally the indices of the f’s can range from 1 to infinity. We suppose 
that the neurons of each type are spread at random through the net. Then 
we shall have the spread of excitation in Jt governed by the following 
equation 


aD Alt expe Evi ()]-« (32) 


However, the following identities hold: 
1— e-“F, (ay) =1-—e% 
1—- e VE, (ay) =1— en Ces 


(33) 


ay ay ay 
1— amt) | (ay) 2 i= @ Set ce, (Paribas oo rent yt 


Cars 


Moreover, if / is the highest threshold in N, we have 


Spa eah (34) 
eit 


Multiplying the right sides of (33) by appropriate factors f; and adding, 
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we obtain 

d = ay? —ay 4 Canyons aa 
aie F, | é ay — Fry 1 é tS peed ee (hee ie y, (35) 
where 


R= pet - (36) 


The steady states are obtained as before by setting dx/dt = 0. 

We shall investigate the special case where two classes of neurons are 
involved with thresholds # = 1 and h# = 2 occurring with frequencies 
fand 1 — f respectively. This situation may be interpreted as involving a 
net with all neuron thresholds equal to 2 but with a fraction f of the neu- 
rons subliminally stimulated. The steady state equation is 

1—e-7— (1—f) aye7—-y=0. (37) 
We have seen that for f = 1, a > 1, equation (37) determines exactly one 
positive stable steady state. For f = 0, a must exceed a certain value for 
a positive stable state to exist. Suppose the condition is satisfied. Then the 
threshold y; and the steady state fraction y2 are continuous functions 
of f, given implicitly by equation (37). By varying f continuously, the 
threshold yy; and the steady state v2 may be made to vary continuously 
also. Physiologically, this may mean that the activity of a center C, 
whose axones innervate the neurons of 9t subliminally, determines the 
overall threshold (“ignition point”) of Jt and the steady state excitation 
which 9 will attain if stimulated sufficiently by an initial outside stimulus. 

Such a mechanism may possibly underlie a “‘preset”’ response, where 
an “‘estimate’’ of the desired “‘intensity’’ of some act seems to be made in 
advance, so that the act can be initiated by a trigger stimulus and then be 
performed automatically with the desired intensity. Possible examples are 
estimating the distance to be jumped, the pitch of a tone to be sung, etc. 
For experimental results with a possible bearing on preset responses, see 
the references given in the paper mentioned above of R. Mayne. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Laboratories 
and the University of Chicago. 
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A theory is presented for the study of diffusion in heterogeneous tissue-like 
structures. It is applicable to a common type of measurement in which the 
change of the amount of substance remaining in the tissue is determined as the 
substance diffuses from the tissue into an adjacent medium, for instance, Ringer’s 
solution. The main objective of this paper is to obtain a method for the calcu- 
lation of the diffusion coefficient in the intercellular space and of the permeability 
coefficients between this space and the cells, based on the type of measurement 
mentioned above. Although the fundamental ideas upon which the theory is 
based are applicable to any type of tissue, the formulae derived are limited to the 
case in which the cells form a flat bundle of parallel fibers. The theory is applied 
to the experimental] results of E. J. Harris and G. P. Burn on diffusion of sodium 
in the sartorius muscle of the frog. We find that if we know the ratio of the 
cellular and intercellular volumes of the muscle the ratio of the equilibrium con- 
centrations of sodium outside and inside the cells can be determined. A very 
simple mathematical analysis of the experimental relation between the amount 
of substance diffusing out of the muscle and the time of diffusion gives us this 
ratio. The ratio of the equilibrium sodium concentrations in the case of the sar- 
torius frog muscle is between about 10 and 30, depending on the muscle used. 
The same mathematical analysis makes it possible to obtain the permeability 
coefficients of muscle fibers through simple calculations, if their sizes are known. 
The permeability coefficients for the experimental] work mentioned above using 
sodium are 1.25 to 11.5 X 10-8 cm/sec for the flow into the fibers and 3.2 to 
16 X 10-7 cm/sec for the flow in the opposite direction. The determination of 
the diffusion coefficient in the intercellular space is more laborious and yields 
only an order of magnitude: 10~-® cm?/sec. 


The objective of the paper and its fundamental equations. Harris and 
Burn have published an experimental and theoretical study on diffusion 
in tissue in which an attempt is made to characterize the properties of a 
tissue by considering the process of diffusion in the intercellular space 
separately from the permeability through the cell membrane. The experi- 
mental technique used by these authors (1949) consisted of following the 
time dependence of the total amount of a substance existing in a tissue 
and diffusing outside of it into a solution of known concentration. The 
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diffusion coefficient in the intercellular space and the permeability coefh- 
cients between the latter and the cells were obtained from a mathematical 
analysis of that time dependence. This seems to be the first time that in- 
stead of a determination of an average diffusion or permeability coefficient 
of the tissue as a whole a more accurate analysis reflecting its hetero- 
geneous structure has been attempted. 

It is the purpose of the present work to develop a theory based on the 
fundamental laws of diffusion and permeability. We shall also take into 
account the heterogeneous structure of the tissue and show the applica- 
bility of the theory to the type of experimental result obtained by Harris 
and Burn. 

Consider an idealized geometry of a tissue slice (Fig. 1) of unit thickness 
in the direction perpendicular to the diagram. Let the length of the slice 
in the plane of the diagram be infinite, let the cells all be of equal size and 
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Ficure 1. An idealized geometry of a tissue slice upon which the theory is based 
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shape and have a fiber-like prismatic structure of cross-section cccc with 
axes parallel to each other arranged in a regular fashion in parallel rows 1, 
2,....,%,...., 2N. Consider the case in which the diffusing substance 
is initially distributed in a homogeneous fashion throughout the whole in- 
tercellular volume. Let the concentration of the substance be the same 
initially in all cells. At ¢ = 0 the tissue is immersed in a solution and the 
substance diffuses into the surrounding medium through the two end 
layers nm = 1 and m = 2N. We assume that the total number of layers of 
the tissue is even. Then the study can be limited to one-half of the tissue, 
since an axis between the layers N and N + 1 is an axis of symmetry 
(Fig. 1) through which no diffusion occurs. Consequently the diffusion in 
one-half of the tissue (Fig. 1), ie., for 7 = 1 to m = N, is the same as 
would occur if this part of the tissue had its layer m = 1 in contact with 
the solution and the layer x = N in contact with a solid into which there 
could be no diffusion. Thus we see that the determination of diffusion in 
one-half of the layer shown in Figure 1 actually solves two different prob- 
lems, both corresponding to the previously mentioned measurements of 
Harris and Burn. The theory developed here could also be extended to 
the case in which the total number of layers of the tissue, as represented 
in Figure 1, is odd, if we consider that the substance flows from the middle 
layer into both adjacent layers at an equal rate. 

The idealization of the tissue as represented in Figure 1 makes the struc- 
ture a periodic one, the strip between the axes aa representing a full period 
and the axis bb being an axis of symmetry within the period. Naturally 
the ‘‘cells” of Figure 1 are rectangular only for convenience of discussion ; 
mathematical considerations do not imply any particular shape for the 
cell. The axes aa are also axes of symmetry; consequently, no substance is 
diffusing across the lines aa and bb. The diffusion goes on in each strip 
aabb independently of the other strips. In general, the diffusing substance 
is distributed throughout the whole tissue. We assume that insofar as the 
overall diffusion process is concerned the distribution of the substance in 
the intercellular spaces may be approximated with sufficient accuracy 
by using the average distribution over the intercellular space associated 
with any cell instead of the actual distribution. This assumption can be 
formulated in the following manner. 

Consider a cell cccc of the mth row. Let C, be the average concentration 
in one-half of this cell, which is determined by the boundary of the latter 
and by the axis aa. Because of symmetry C,, is also the average concentra- 
tion in the whole cell. We assume that the behavior of the cell in the tissue 
is described with sufficient accuracy by the magnitude of this average con- 
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centration alone. We consider each cell to be surrounded by an adjacent 
intercellular space with its external boundary at inl (Fig: 1). LetiT, be 
the average concentration in such a space of the wth row. The average 
concentration in each half of the intercellular space bounded by aa is 
also I, for reasons of symmetry. Whenever concentrations have to be as- 
signed to well defined points, we will imagine the concentrations I, exist 
at certain points P, inside the intercellular spaces. For reasons of sym- 
metry there are two such points around each cell. On the basis of this 
scheme diffusion flow connects each half of a cell directly only with the 
adjacent half of the intercellular space. Consequently the change of the 
amount of substance in the cell is, per unit of time, 

AdC,, 


t= S(hiIn— hn), = 1,2,....,N, (1) 


where A is the cross-sectional area of the cell, i.e., the volume of the cell, 
since the tissue is assumed to be of unit length in the direction perpendicu- 
lar to the diagram; S is the area of the cell membrane, &; and k, its permea- 
bility coefficients. The subscript 7 is limited to 1,...., NW because we 
consider only one-half of the tissue up to the axis of symmetry. We shal] 
now consider that part of an intercellular space of the mth row which lies 
between the axes aa and bb. Three diffusion flows exist between this space 
and its surroundings as indicated by arrows in Figure 1. Two of these 
flows are from or into the adjacent intercellular spaces of the rows n — 1 
and ” + 1, one is from or into the surrounded cell. The average concen- 
tration gradient in the diffusion from the intercellular space surrounding 
P,1 to the intercellular space surrounding P, can be taken with good 
approximation to be 
ee, W=e2,....,V—1), 
where L is a length which does not differ by very much from the length 
of an arc like P,P, (Fig. 1). This is about a linear dimension of the cell 
plus the intercellular distance. We assume that L is the same for any two 
adjacent rows, that is, L is independent of m. The number of molecules 
which diffuse per unit time from the part of the intercellular space sur- 
rounding P,_; to the part surrounding P, can be written in the form: 


Tuer ey nh 
DW 
where W is an average width of this diffusion flow and D is the diffusion 
coefficient in the intercellular space. The length 2W does not differ greatly 
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from the average distance of two neighboring cells. Using these symbols 
the change of the number of molecules in the intercellular space bounded 
by iii and cccc is, per unit of time, 

ae = aa mento (Rein — Ria) (N= Diy we IN —4) (2) 
where B is the cross-sectional area of that intercellular space, 1.e., its vol- 
ume, because we consider a slice of tissue of unit thickness. In numerical 
calculations B is often taken to be somewhat smaller than that volume in 
order to take into account the solid substances that it contains into which 
practically no diffusion flow occurs. Since no substance is flowing across 
the axis of symmetry we have for m = N, instead of equation (2), 


Dien ee Ue eee (3) 


We assume that the concentration of the diffusing substance is negligible 
in the external solution. Then the concentration gradient between layer 1 
and the boundary of the tissue is /;/(L/2) and, instead of equation (2), 
we obtain for m = 1: 


Ben OY 31) —S(krli— ki). (4) 


The diffusion problem outlined corresponds to the following initial condi- 
tions for J,(¢) and C,(t): 


TRO =e, (S$) Ose. kG, (0) —k71, (0) | forn=1,..2.4N. (©) 
t=0 


These conditions imply an initially constant concentration J in the inter- 
cellular spaces and a constant concentration in the cells in equilibrium 
with those spaces. The cellular concentration C, in any of the layers can 
be calculated as soon as the intercellular concentration J, in the same 


layer is known. In fact, from equation (1), we have forn = 1,...,N: 
Edm aft t) dl 6 
OGM =EIth sf eOhOd, (6) 


where e(#) = exp (k.St/A). A 

Take a portion of the tissue of wit thickness consisting of NV layers 

(n= 1 Norn =n+1,..., 2N) with M cells and M intercellular 
ma es 


spaces in each layer. For this portion of the tissue let 


a = MNA = total cross-sectional area of the cells = total cellular 
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volume, 


8 = MNB = total cross-sectional area of the intercellular spaces = to- 
tal intercellular volume involved in the diffusion process, 


« = MNS = total surface of the cell membranes, 


\ = NL = average length of the diffusion path within the intercellular 
spaces across the whole portion of the tissue from one boundary to the 
axis of symmetry. The quantity 2) is likely to be somewhat larger than 
the size of the whole tissue in the direction aa. 


X, = aC,/N = amount of substance in all the cells of the mth row, 


Y, = BI,/N = amount of substance in all the intercellular spaces of 
the nth row, 


n=N 
OS Se (X,+ Y,) = total amount of the substance in the portion of 


n=1 


tissue considered. 


Using these symbols equations (1) and (2) can be transformed into: 


dl (005) ig ON, pe (7) =e Vis 
ave (7) 
di = (VY,1—-2V,+ Var) —eY +eiXa; (1 =2,....,N—1), 


where ¢; = k.a/a, go = kyo/B, 9 = 2DWM/(Bd), — = oN, X,(0) = 
(v2/¢:)V/N, Y,(0) = Y/N, Y = BI is the total amount of substance 
existing initially in all the intercellular spaces, and (¢:/¢:)Y is the total 
amount of substance existing initially in all the cells of the portion of tis- 
suen=1,...,Norn=N-+1,...,2N. Equations (3) and (4) yield 
two other relations obtainable from equation (7) if we put Vo = —Vj, 
Vysi = Vy by definition. 

2. Solution of the equations; total amount of the substance in the tissue. 
To solve the system of equations derived at the end of the previous sec- 
tion we shall apply the Laplace transformation, that is, associate with 
each function F(é) its Laplace transform f(s) as defined by 


f(s) =f oP) i. 


We will use the following type of symbolism: 
f(s) =QiFO}, FG) =2 {Ff (s)}. 
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If we put #,(s) = U{X,()} and y,(s) = &{Y,,(é)} we obtain: 


S%n — Xn (0) = GoVn — G1%n (GUS 2h ce, IV) 
SVn — Ve (0) = E (Vn—1 ay 2 n+ Vn-+1) (8) 
—G29n F G1%p ee Oe eet ie 
The last equation is valid for 7 = 1 and m = N if we define yp = —y, 
and yv41 = yy. Taking into account the expressions of X,(0) and Y,,(0) 
we obtain for n = 1, 2,..., WN, using the above definitions of yo and 
INH) 
X, = 22 fon a 1 — N-1 
21 s+ [ + oe (5s + 41) ~ EWS Ve N Wo) (9) 
= E (nat — 2In 1b Yn) 
The solution of equation (9) form = 1, 2,..., N is: 
a V As ont g2Nt1—n 
n= Syl 1-? Gender | OY 
where zg = u + Vv? — 1; or 
i 
s+3). ole Usha) ae 


The quantities « and y, do not change if z is changed in them into 1/2, 
because the coefficients of y,1 and y,4: in the difference equation for y, 
are the same. 

Let Q(é) be the amount of substance existing at any time ¢ in the por- 
tion of tissue from ~ = 1 to m = N. The Laplace transform of Q(#) is 


g(s) =210(} Si +») = = (142) (1-205) 


RO) =(14+2) Ute(sto Ir 


2 PNG GG) 


d=. (31) aE); 


(12) 


The total amount of substance initially existing in the portion of tissue 
considered is 


Q(0) = y (1+). 


Consequently G(¢t) = Q(4)/Q(O) is the fraction of the initial amount of the 
substance which at the moment t still exists in the tissue. We obtain for the 
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Laplace transform of G(¢): 


1—h(s) 
neoteeor == 
An easy check of the calculations made up to this moment is possible. The 
function u tends to 1 as s > 0; consequently z tends also to 1. Therefore, 
h(0) = 1 and, we find, by an application of a Tauberian theorem of the 
Laplace transformation (Doetsch, 1937, p. 208, theorem 2)2 


oe {ats ty when to. 


Consequently G() = 0, as it should because at ¢ = ~ no substance is 
left in the tissue. When s— ©, uo and z— . Consequently 
h(s) +0 and g(s) behaves at infinity as 1/s. Therefore G(O) = 1, as it 
should by the meaning of the function G({) itself. 

We now calculate ¢"'{7(s)}. The poles of T considered as a function 
of z are the zeros of the denominator of T that are not zeros of its numera- 
tor. Consequently, these poles are roots of the equation 


ge Lise OG; Cis) 
i.e., they are the following values of 2 = z,: 
Dip 
2, = exp (10,) where 0, = ( eas With 1 = 1, 25 43 Vee 


In order to discuss the singularities of T considered as a function of u we 
need the following theorem: If a function F(z) has a simple pole at z = zp, 
and if w = uw, is a simple root of the equation 2(w) = zn, where z(u) is an 
analytic function of u at u = u,, then u = u, is also a simple pole of 
F{z(u)] considered as a function of w. The proof of this theorem is easy. 


In fact, 
K 


Ey ee 


+6(2z), 
Z— 2, = (U— Uy) V (u) , 


where ®(z) and W(w) are analytic functions at z = z, and u = wu, respec- 


tively, and K and W(w,) are two finite constants different from zero. 
Therefore 
ERG 


(u—u,) Fl 2 (uw) ] Saey 


+ (u — Up) & [Lz (u) ]. 
Since P[z(w,)] is finite, W(w,) and K are both not equal either to zero or 
infinity, we see from the last equation that wu = u, is a simple pole of 
F|z(u)] with the residue K/V(u,). This proves the theorem. 

To each value of g = g, corresponds by the first equation of (11) one, 
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and only one, value of u which is u, = cos 6,. The function 2(u) is analytic 
at all ~ = u,. Consequently u = u, are simple poles of T considered as a 
function of u. Since 

Cos 8,, = COS Aon—n+1 5 (15) 


the poles of T considered as a function of z [see equation (14)] can be 
grouped in pairs so that each pair gives rise to a single pole of T con- 
sidered as a function of w. Consequently, in dealing with this latter func- 
tion we can limit the range of the subscript n ton =1,..., N. It may be 
easily seen that beyond wu = wu, there are no other singularities of T since 
these could arise only out of the singularities of the function z(w) which 
are u = +o, +1 and correspond to z= 0, +~, +1. A direct inspec- 
tion shows, however, that these are not singularities of 7. 

We proceed now to discuss T as a function of s. The relation between s 
and wu is given by the second equation of (11) which can also be written 
in the form: 


+ (git getr&) strmée=0, (16) 


where v = 2 — 2u. Two values of s which we call s, correspond to each 
value of u = u,. They are the roots of the following equation in s,: 


Sat (eit gotiné) sn tmeré =0 ; (17) 
where v, = 2— 2 cos 6, andu = 1,..., N. The roots of equation (17) are 
always simple, both real and negative. The larger of these roots will be 
indicated by s}, the smaller by s;. It is seen from the second equation of 
(11) that the function u(s) is analytic at all s = s,, and, therefore, 
s=s} and s=s, are simple poles of T considered as a function of s. 
It may be easily seen that beyond s = sj, s; there are no other singulari- 
ties of T, since these could arise only out of the singularities of the func- 
tion u(s) which are s = —g1, + © and correspond to wu = + ~. But the 
latter are not singularities of T. 

We can now apply the Heaviside partial fraction expansion as follows: 


eS) > Se, 


Sch = 
See Sis ante 1 


2 in il 5S — Sn 
ees uP exp (5,1) Bess ares| PN: 
Since 22” = —1 and, consequently, "1 = —2,', we obtain: 


2 on 
f,= WN? z <> ras Oa (sal) (FZ mh 
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If we put p=gits, equation (5) becomes: 
pi (g1—- go. —vE) P-—Gige=9 : 


From this we find: 
ds dp dv du _ byl |! 
az dy du dz 2 $102 


From the last expression of 7,, we obtain: 


sik 
T, =20(1+ 9%) exp (Syl) - (18) 
From the expression of &"{7(s)} and of g(s) we obtain, by a routine 
application of the procedures of the Laplace transformation: 


G(t) =1 


n=N 


oy. 29 pee go -exp( 5,0) a 
$3 git go Pn : 


’,=ti.t, n=l 14 218 
£, 


where pf = g1 + s; and p, = ¢1 +5; are the roots of the following 
equation in py: 


Pn — (G1 — G2 —mE) Pn — G1Ge =0. (19) 


From equation (19) we see that p; > 0, p, < 0. Consequently, if we put 


Ri = V (¢1— ¢2 —mE)? + 4ei@: , 


we easily obtain for pa = pi, pr: 


Using this relation it may be seen that the term in exp (—¢ié) in G(d) 
vanishes. Consequently, if we put 


et = exp (s*?), 6, = exp (s7i), 


we can write: 


G = 14 S(t S)a-o -(4+#)a-«) 


The condition G(0) = 1 can be easily checked for this relation. The rela- 
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tion G() = 0 gives the following identity: 
n=N 


2 On 
>, cosec? = 2N?. (20) 


n=1 


If this is taken into account the expression of G(t) can be further sim- 
plified: 


‘ 2 n=N : 
G (1) a Pi Tn + 20n (21) 


it go net R, f 
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Ficure 2: Analysis of diffusion of sodium out of a frog sartorius muscle (muscle W 46). 
Fraction of substance remaining in the tissue against duration of diffusion. 


where 
Dee ap gee ees 
n s7 st n? n Ca gt 
n n nr n 


This is an expression for the fraction of the initial amount of the substance 
which at time t exists im the tissue. 

3. Determination of diffusion and permeability coefficients from experi- 
mental data. Figure 2 represents the results of an analysis of a set of data 
kindly supplied by Drs. Harris and Burn. The circled points joined by a 
fully drawn curve and one surrounded by a square are experimental 
values on the diffusion of radioactive sodium from a frog sartorius muscle 
immersed in Ringer’s solution. The thickness of the muscle was 800y; its 
mass, 75 mg. The ordinate represents the fraction of radioactive sodium 
remaining in the muscle after immersion plotted on a logarithmic scale. 
The abscissa gives time in minutes. This set of data was analyzed with . 
the objective of determining the diffusion coefficient D of sodium in the 
intercellular spaces and the permeability coefficients k; and k, of the cell 
membrane. It is apparent from the diagram that the function G(¢) tends 
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asymptotically to be exponential and this fact has been found very con- 
venient for the calculation of the constants ¢: and ¢2 which determine the 
permeability coefficients. Since the expression of G(/) is of the type 
DA, exp (s,t) with A, independent of ¢, the simple exponential shape 
of G(x) for large ¢ is due to the fact that from a certain ¢ on the term con- 
taining the largest s, is much larger than the sum of all the remaining 
terms. To analyze this fact it is necessary to examine the location of the 
quantities s, or p, which are the roots of equations (17) and (19) respec- 
tively. From equation (19) we see that p; and p, have opposite signs; 
consequently, pt > 0, p, <0. Therefore; > —qi, 5; < —¢u. Since 
st has the same sign as s, we can see, using equation (17), that both s; 
and s, are negative and can be represented diagrammatically as in Fig- 
ure 3. Differentiating equation (19) with respect to 2 and taking into ac- 


=o «<—— S$, S, Ss, $s, -90,<— Sa SSeS o - 
See = = 2 + tte ae 
FIGURE 3 


count the expression of v, we obtain: 


Haale a Ep, sin @, . 
Consequently d| py |/dn < 0 and d|p,|/dn > 0. Therefore, as ” increases 
Pn approaches 0 and pf, moves away from zero, as indicated in Figure 3. 
Since the largest m is N and #y > 7 as N — = it can be easily seen 
from equation (19) that py > 0, py ~ —@. 

For large ¢ the terms in G(¢) containing exp (s;#) are negligible with 
respect to those containing exp (sft). A preliminary analysis of the data 
of Figure 2 showed that 2 was of the order of magnitude of 10-5 sec 
whereas p,; was between 10-* and 10-8 sec. Consequently, as a first 
approximation, p; could be neglected with respect to g, and we put 
a ~ Q in the expression of G(). In addition to this, the same analysis 
showed that gig, was of the order of magnitude of 10-° sec-? whereas 
Pné — (~1 — ge) was of the order of magnitude of 10-° or 10-3 sec—. 
Consequently 


R, = |mE— (gi — g2) |. 


But, whereas »,£ appeared to be of the order of magnitude of 10 or 


DIFFUSION IN MUSCLE 57 
10~* sec, 1 — yg: = ¢ was of the order of 10-4. Therefore: 
5 ey hs 
R, ~mE = 4N2o@ sin? 3° 


Since p; was of the order of 10-* to 10-8 sec— we put sf ~ —g,. In this 
way we obtained from equations (21) and (20), for large ¢, 


we eT (22) 


This procedure has been found quite accurate in reference to the data of 
Figure 2. The curve between 130 and 240 min is practically a straight line 
which gives us, using formula (22), 


gi = 1.426 X 10~4 sec! 


’ 


go= 3.124% 10 sec! . 


At ¢ = 240 min the curve in Figure 2 seems to pass into another straight 
line which has been interpreted as a consequence of some change occurring 
in the muscle after the first four hours. This is analyzed in section 7. For 
the present the data of Figure 2 above ¢ = 240 min may be disregarded. 

The determination of ¢ is somewhat more laborious and requires an 
actual fitting of the curve G(¢) with its theoretical expression. A crude ap- 
proximation to the value of g can be obtained by comparing the theoretical 
expression of the slope of G(¢) at ¢ = 0, which is 


pee (23) 
git oe 


with the one obtainable from the experimental data. It will be shown 
shortly that 5 is a reasonable value to assign to N for the muscle in 
question. Estimating by numerical differentiation the slope of the graph 
at t = 0 in Figure 2 as —4.54 X 10-® sec~ and using the previously 
determined values of g; and ¢, formula (23) gives ¢ ~ 5.7 X 10-* sec". 
This has been found quite good but somewhat too small to fit the whole 
curve. A better value of ¢ has been obtained by using the exact expression 
of G(#) and determining by trial and error that value of y which gave the 
experimental value of G at ¢ = 25 min. Two or three trials accompanied 
by an interpolation gave us: 


pe (7a 10 “sec. 


The plus signs in Figure 2 and the point surrounded by a square represent 
the values calculated with the above ¢ by means of the exact expression 
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(21) of G(t). Beyond? = 50 min the theoretical curve smoothly approaches 
the experimental one. 

Although for reasons of mathematical simplicity Figure 1 has been used 
to derive the equations of the theory, Figure 4 represents somewhat more 
closely the arrangement of the fibers in the muscle considered. Three fibers 
of circular cross-sections placed in two adjacent layers are represented in 
this figure. The centers of cross-sections of the fibers form an equilateral 
triangle. The minimum distance between two neighboring cells 2KW is 
somewhat smaller than the width 2W of a single intercellular diffusion 


Ficure 4. Schema of fiber arrangement in the muscle 


flow, i.e., K < 1. The following relation holds between the radius r of the 
fiber, the thickness / of the tissue, and the number 2N of layers: 


2N(r+KW) V3=1. (24) 


If we denote by » the fraction of the volume of the tissue occupied by 
the fibers, i.e., the ratio of the areas of three circular segments like ABC 
to the area of the triangle ADE, we obtain: 


rre=2V73n(r+KW)?. 


Boyle et al. (1941) give for » a range from 0.87 to about 0.91. Mayeda 
(1890) gives for y an average value of 42u. On the basis of these data the 
number 2N of layers in a tissue of J = 800 can be calculated, since the 
previous two relations give: 
2 P 
2N)?=—e 
GN arvV3 2° 


From this we find 2N = 10.77 for n = 0.87. We take 2N = 10 which cor- 
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responds to r = 45u with the same value of 7. This is still a reasonable 
magnitude for an average size of the fiber. The distribution of the radii of 
the sartorius muscle in frog shows, according to Mayeda (loc. cit.), a 
sharp peak around a mean of 42y, but the full range goes from 10 to 76y. 
The above value of r implies, using formula (24), KW = 1. Accepting 
the values of ¢, and ¢ previously obtained and putting L ~ ar, S ~ 2ar, 
we obtain the coefficient of permeability: 


how Or = 3.2 X 10-7 cm/sec, 


TABLE I 
Muscle 5 Thick- 
i Weight 1 2 ke kr 

ao Mg. ns Sec.7} Sec.-1 Cm. Sec.-1 Cm. Sec.71 Ke/kr 

N 46 iD 0.08 1.43107! 3./2X< 105° SAS e 1 26>S10Re e255 

N 67 ts: 0.08 3230<10m4 IO xLOm4 7.A3X10-* 5.91 X10 |) 12°56 

N 69 90 0.1 S252 1 0e SHO Ons 2 Oe 1 Oe ea StS LOSS 1 1055 

N 70 90 Oeil diel 10m Dro0X<tOwen eto. Ome<On 4 9.44x10-8 | 17.0 

Ni71 85 0.08 3.8710 8.88X10% Sail Oma 2.99X10-8 | 29.2 

N 72 85 0.08 AAAS 10s 3.04X10~4 OF25 10m 11052251052 9.06 


All experiments done at 18° C. 


which is about three times smaller than the values indicated by Harris and 
Burn, and 


ie easel = 
ky= gor TAL SS Lao S< 1G) 8 cm/sec 


which is about 10 times smaller than a value of Harris and Burn. In this 
way the ratio of sodium concentration outside of the cells, which is in 
equilibrium with the concentration inside, is ~ 26 for the data used against 
a range of about 3.5 to 7 indicated by Harris and Burn. For the diffusion 
coefficient in the intercellular space we obtain: 


D =0.8 (xN r)? (qn — 1) eo K X 10-5 cm?2/sec. 


The coefficient 0.8 in this formula takes into account the solid substances 
existing in the intercellular spaces into which no sodium diffuses, K is the 
ratio of the minimum distance between two neighboring cells and the 
average width of the diffusion flow between these cells. Consequently 
K <1andD can be estimated as being about the order of magnitude of 
10-® cm?/sec, which is in agreement with the calculations of Harris and 


Burn. . 
Table I gives results of calculations of the permeability coefficients 
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of a few other muscles, the data for which have also been kindly sup- 
plied by Drs. Harris and Burn. The muscle NV 46 is the one whose analysis 
was just described in detail. 

4. Amount of substance inside and outside the cells. If Q,(t) is the amount 
of substance existing at t in all the intercellular spaces, we have: 


= SY. 0. 


The same amount expressed as a fraction of the amount initially existing 
in all the intercellular spaces is: 


For the Laplace transform of the above we have: 
beef 
g1(s) = {Gr } =——. 


Applying the same procedure as the one for g(s) we obtain: 
n=N 
Tn 
Gr(V) =o DR (25) 


Similarly for Q.(/), the amount of substance existing at t in all the cells, 
we find 


WO = Sek 


If we express this as a fraction.of the amount initially existing in all the 
cells we have: 


n=N 


6. =o ate. 


For the Laplace transform of this we have: 


g.(s) =21G, tt) } =e 
1 


From here we obtain: 


n=N 


On 
G, (t) =2¢¢; Re (26) 


The following relation can be easily checked with the above equations: 


%G 14 (2 
Cit al 1+(2)a. (27) 
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5. Distribution of the substance among the various layers. The amount of 
substance existing in single layers can be calculated without additional 
difficulty. From y,,/(Y/N) = [1 — h,,(s)]/s, using 
gmt g2N+1—m 


(TF 2) (TF 29 ’ 


one can obtain the amount of substance existing in the intercellular Spaces 
of the mth layer. Since 


hy (s) =2 


n=N 


L71 {hm (5) } =2NG >” | fn| A mon EXD (Sn) 
n=1 


where 
{ _ cos (m — 1) 6,,—cos m6, 
eae R, ‘ 
we have 
n=N 
Va e= 2s A winte (28) 
n=1 


The relation Y,,(~) = 0 gives us the identity: 


n=N 


>» sin (m Gee 3) 6, 
A : 
n=1 Ss DEES oa 
sin 5 N 


The amount of substance im the cells of the mth layer can be calculated 
from the following equation which is easily obtainable from the relation 
between C,(t) and J,(¢): 


: t 
Xm () exp (gl) =2 4 os [Vn (0 exp (ord) dl. 
0 


Pi N 
We have: 
n=N 
X mt) = 2e—2V DY) A mindn- (29) 
n=1 


Since for large ¢ 
pac hee 202) 


WT, 0; oA a 


as the diffusion goes on the amount of substance in the intercellular spaces 
becomes negligible with respect to that remaining inside the cells. Taking 
into account also that 


R,~ 4N’%¢@ sin? - 


it is easy to show that the amount of substance remaining in the cells tends 
to be approximately the same in each layer as ¢ becomes large. 
6. Formulae at the beginning of the diffusion process. For small i we ob- 
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tain by power series expansion of €,: 
KS 


= —R,t A att ae 
™ NY @ 4 
R R 
a Ny 
Ph Ga 


Therefore, if we take into account the trigonometric identities previously 
derived, we obtain for the amount of substance existing im the whole tissue 
at the time t, divided by the initial amount: 


GU) a ee (30) 
git ¢e 
This amount is distributed as follows among the intercellular spaces and 
the cells: 
Geen ae (31) 


G. (t) =1— Neel’. (32) 


I 


Comparing these three relations with each other and keeping equation 
(27) in mind, it is seen that g, and ¢, appear in the first order term of G(¢) 
only because at ¢ = 0 the distribution of the substance among the inter- 
cellular spaces and the cells depends on g; and ¢g». At the very beginning 
the rate of diffusion out of the tissue does not depend on the permeability 
coefficients of the cells. 

7. Conditions at the boundary between the tissue and the solution. An in- 
spection of Figure 2 shows that the theoretical fit is not as good as one 
would desire for small ¢. It cannot be improved, however, since it would 
require a lowering of ¢ which would make the fit worse for larger ¢. The 
question arises whether the fit could be improved by taking into account 
the fact that the concentration of the substance at the boundary may not 
be so small as to be taken equal to zero, particularly at the beginning of 
the diffusion process. The objective of the present section is to analyze a 
possible effect of this situation. 

Let us consider in the solution two layers adjacent to the tissue, and 
let us number them » =.0 and m = 2N + 1 (see Fig. 1). Let us assume 
that the concentration of the substance in the solution can be made equal 
to zero at points like Py within these layers not very distant from the 
boundary of the tissue. Consider a part of the layer n = 0 between two 
axes aa and bb and the diffusion flow between it and the adjacent inter- 
cellular space of the layer » = 1. We consider Wy to be an average width 
of this diffusion flow (W» is somewhat larger than W, but smaller than V), 
Dy to be the average diffusion coefficient of the substance considered in 
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the intercellular space and in the external solution. Then the number of 
molecules which diffuse per unit time within aabb from the intercellular 
space under consideration to the adjacent part of the solution (Fig. i) is 
approximately : 
D.Wol, HDWI; 
jin, a a 


where Ly is the average length of this diffusion flow, that is, about the 
length of an arc like P,P. The coefficient H defined by the above equa- 
tion is introduced for the sake of mathematical simplicity. Consequently, 
we find: 

Bdl,_2DW 2H DW 


ae ae i I,—1,) — 7 


dt IE 4 tp (ep 8 4) 


Since it is likely that Do, Wo, Lo are of the same order of magnitude as 
D, W, and L respectively, we will put H = 1 tentatively, which reduces 
the above differential equation to equation (2) taken for » = 1 with 
I, = 0. With this definition of J) and with the additional condition Jy; = 
Iy expressing the symmetry of the tissue with respect to a middle axis in 
the tissue, equation (2) can now be used for m = 1,..., N, and equa- 
tions (3) and (4) can be disregarded. Equation (2) is solved in a similar 
manner as before and g(s) still has the same expression with the only 
difference that 
Find Ot ele ee eT gee eal eaten 

WN { 21) (2783 1)? : 2N+1 ° 
instead of equations (12) and (14). Equations (16) and (17) remain un- 
changed and the right-hand side of equation (18) has a factor y cos? (6,/2) 
instead of 2v, where y = 4¢N/(2N + 1). The same replacement of 29 
has to be made in equations (21), (25), and (26). Equation (22) is un- 
changed, so that this theory yields the same permeability coefficients ky 
and k, as the previous one. The theoretical slope of the curve G(é) at 
t = 0 is now only 3 of the previous one. Equations (28) and (29) are now 
changed by the replacement of 29 by y and by a new meaning of Ann: 


__ sin 6, sin (m 6,,) 
msn J ie 


The trigonometric identities which hold for the present angles 0, are: 


n=N rs) 
Dcot®s= AN+1)N, 
n=1 
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instead of equation (20), and 


n=N 6 
i eae ie 
BPE MGA 2 N+ 


The coefficients of ¢ and # in equations (30) to (32) are now only one-half 
as large as before. 

A fit of the experimental data of Figure 2 by means of the present theory 
has been done in the same way as before. The points marked by multipli- 
cation signs in Figure 2 are points of the theoretical curve which goes 
through the experimental point at ¢ = 25 min, has the same values of 
v1 and ~ as before and ¢ = 7.85 X 10-4 sec. It is seen that although 
the new fit is somewhat better than the previous one for larger ?’s, it is 
somewhat worse for smaller ¢’s. The difference between the two values of ¢ 
is about 5%, which is completely irrelevant in view of the fact that ¢ is 
used only for the determination of the diffusion coefficient in the inter- 
cellular space, and this determination is affected by other uncertainties 
anyway. An interesting fact is that both theories lead to exactly the same 
permeability coefficients of the cell. Thus it does not seem possible to ex- 
plain in this way the slight divergence of the theory from the experi- 
mental data of Figure 2, particularly at small ¢. However, the difference 
between the experimental data and the theory for small ¢ may be partially 
accounted for by diffusion from the muscle in several ways not taken into 
account by the theory. We have considered the muscle to be an infinite 
sheet, whereas it actually has a small roughly triangular area instead. Thus 
diffusion actually occurs through all the edges of the muscle. Since approxi- 
mately 10% of the total volume of the muscle lay within edges at a dis- 
tance equal to twice the thickness of the muscle, the more rapid diffusion 
measured than given by theory is not surprising. 

The dotted curve of Figure 2 represents a fit of the experimental data 
on muscle N 46 made by using the data above, ¢ = 240 min, for the deter- 
mination of g: and g (¢g: = 1.12 X 10-4 sec, g = 1.86 X 10-5 sec). 
The curve, calculated on the basis of the theory of this section, has been 
made to pass through the experimental point at t = 25 min which required 
g = 6.16 X 10~* sec“. It is seen that although the fit is very good above 
¢ = 130 min, the divergence is very strong for lower ?’s. As previously 
mentioned, this is the reason why the data above ¢ = 240 min have been 
disregarded in final analysis. 

This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of The University of Chicago. 
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The outputs of nervous systems (as expressed in motor activity) are viewed 
as mathematical transformations on the inputs which enter via the sensory 
nerves. Simple nerve-ganglion models are exhibited which theoretically account 
for the arithmetic computations necessary to expedite such transformations. 


Introduction. As a further development of the statistical approach to 
the theory of the central nervous system, it would be desirable to develop 
equations that would in some way delineate the complexity and versatility 
of the neural models implied by this orientation. 

Simple nerve ganglion systems, for example, imply a variety of input- 
output problems. The output of such a system is in general a transforma- 
tion of the input which is mediated by the intervening ganglion. It would 
be interesting to know what kinds of transformations are possible with 
ganglia whose only variables are size, distribution of synaptic delays, dis- 
tribution of refractory periods, and distribution of periods of latent addi- 
tion. 

It may be that even such simple systems can exhibit various types of 
abstraction. The output (varying with time) of certain ganglia may, for 
example, give invariant responses to certain classes of time variable 
inputs. 

Arithmetic Computation. If we assume that the function of the brain is to 
perform the proper mathematical operations (transformations) on the 
multi-modal inputs which reach it via the sensory nerves, and if we further 
assume that the arithmetic computations necessary to expedite such 
transformations are. performed by simple neural arrays such as nerve- 
ganglion systems, it is necessary to indicate, at least theoretically, how 
such computation may be performed. 

Addition. We must show, to begin with, that the nervous system is 
capable of responding to the sum of several stimuli of the same modality. 
A “possible” way in which the activity (average firing frequencies) of 
separate axone bundles can be summed is illustrated in Figure 1. 
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Each of the fibers of bundle number 1 enters the ganglion and synapses 
upon only one cell body. The same is true of bundle number 2. Each cell 
body of the ganglion gives off only one axone. These axones emerge from 
the ganglion to form the outgoing bundle (number 3). 

Let f; be the average firing frequency of a typical fiber in the 7th bundle 
(average number of action currents per unit time per axone). Also let F; 
denote the average firing frequency of the whole bundle. 

By definition we may write nf; = F;, where m is the total number of 
fibers in the bundle. The model in Figure 1 is simply a composite of 2 units 
such as that illustrated in Figure 2. 


—_ 


FicureE 1 


FIcure 2 


If we assume that either fiber 1 or fiber 2 is sufficient to fire fiber 3, then 
for low input frequencies fs ~ fi: + fo. 

If we assume that f; and f, are Poisson-distributed then the probabilities 
that fber number 1 or number 2 will fire within a 5 of time (duration of 
refractoriness) are roughly 6f; and df: respectively for a small 6. The 
probability, therefore, that fiber number 3 will fire within the same 
period is given by 

Ofs = dfit+ OF oO iad ts (1) 
and therefore 


$= Jit Je — 6 haa. (2) 


But since nf; = F; it follows from equation (2) that 


6 
= hth o Fils (3) 
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Equation (3) tells us that F’; is somewhat less than the sum of the inputs. 
The per-cent of error E is given by the expression 


6 
= Tihs 
n O fits 
E=—=—— X 100 = —=_ 100. 
Fi +F, Tone < 
Since the maximum error occurs when f; = f. = f, we find 
SES eR Net ®, 


Now if we let 6 = 10~* sec. and f = 200 impulses per second we get 
(max) E = 10% error. However, fibers rarely maintain frequencies of 
200 impulses per sec. for very long. If we let f be somewhat more typical, 
say 25 impulses per second, then (max) E = 1.25% error. 

For normal ranges of f; and f2 we can therefore expect the mechanism 
in Figure 1 to record the sum of F, and F, to within less than 2% error. 

Multiplication. The product of two input frequencies can be recorded 
by essentially the same mechanism as that illustrated in Figure 1. In this 
case, however, the thresholds of the ganglion cells are assumed to be higher 
so that both of the incoming fibers must fire (within the time of latent 
addition o) in order that the outgoing axone should fire. 

By similar considerations we find that roughly 


earn fie (6) 
nN 


The output of such a system is proportional to the product of the inputs. 
Such a model, however, may not be too stable in view of the fact that the 
proportionality constant o/n is so small. This implies that the output, al- 
though proportional to the product of the inputs, is relatively much 
smaller. The behavior of the system, particularly for large 7, may become 
erratic due to random fluctuations in the ganglion. 

A more stable model for recording the product of two frequencies would 
be possible if a ganglion could be constructed such that its output is pro- 
portional to the logarithm of its input. Actually all we would need for 
practical purposes is a ganglion whose output would be approximately _ 
proportional to the logarithm of its input for a substantial range of input 
frequencies. 

Such a ganglion would make possible the construction of the model 
illustrated in Figure 3. 

Subtraction. The mechanism illustrated in Figure 4 is designed to 
record the difference between two stimuli. 
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Bundle A sends axones to nucleus WV and also sends inhibitory branches 
to N>. Bundles N; and N» give off axone bundles C and D respectively 
and also the combined branch bundle E. The activity in A tends to pro- 
duce an activity in NV; proportional to the activity in A. The inhibitory 
branches of B tend to inhibit the activity in V1 proportional to the dif- 
ference between the activities in A and B, so long as (a) is greater than (0). 
When (A) is smaller than (B) the activity in C is zero. Similarly the ac- 
tivity in D is proportional to the difference (B)—(A). 


LnF, 


LnF, +LnF, = Ln(F, Fo) 


FIGURE 4 


Finally, the activity in E is proportional to the absolute value of the 
difference between (A) and (B). 

Ratio. The chapter on Discrimination of Relations in N. Rasheeay S 
book Mathematical Biophysics (1938) shows a simple mechanism which 
records the ratio of two stimuli. Many simple alternatives to this model 
are possible. Only a slight change of nomenclature and argument easily 
transforms Rashevsky’s model to one consistent with the statistical ap- 
proach. 

Differentiation. Models designed to answer more complicated needs 
may require similar or identical information to be “fed” into two or 
more sub-systems at the same time. As has already been indicated in 
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Figure 4, such duplication of information can be achieved by the simple 
branching of the axones in the bundle under consideration. The arrival of 
information at its central destination may be delayed (if some mechanism 
requires such a delay) by having it pass first through one or more ganglia. 
At each such “station” the information will be delayed by about half of 
one millisecond. 

Now with the above it is possible to construct a mechanism which will 
respond to the first derivative with respect to time of the activity in any 
input bundle. This is achieved by first duplicating the information (with 
a time lag as discussed above) and sending the two identical streams of 
information (but out of phase) into the mechanism illustrated in Figure 4. 
The resulting transformation would give the difference between the ac- 
tivity in the bundle at a time / and at a certain constant time earlier. But 
such a difference is roughly proportional to the first time derivative of 
the input as required. 

Repetition of this device could be used for abstracting higher time 
derivatives of the same input. 

Although models such as those discussed in the foregoing are interesting 
as such, it would be much more desirable to derive mathematical ex- 
pressions, ‘‘existence theorems” as it were, which would tell us in a general 
way just how much is possible in the way of models starting with our 
initial assumptions. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Laboratories 
and The University of Chicago. 
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The refractory periods of an aggregate of simple ‘‘counter” neurons are assumed 
distributed according to some probability frequency. The output of the aggregate 
is computed for rectangular and triangular distributions. In particular, itis shown 
that the maximum output of an aggregate with any triangular distribution can- 
not exceed the maximum output of its average neuron by a factor greater than 
2 In 2. This puts an upper bound on the amount of departure from the behavior of 
the average neuron which an aggregate characterized by a certain type of dis- 
tribution can show. 

Next, the aggregate is supposed to be subjected to regularly spaced stimuli. 
Under these conditions, a single neuron will give a discontinuous output curve. 
If, however, the refractory periods are distributed according to some frequency, 
the output curve may be “smoothed out.” A general condition on the distribu- 
tion is derived which makes the output monotone increasing with the input. The 
condition is applied to some special cases. 


If each neuron in an aggregate were characterized by the same parame- 
ters, and if the neurons did not interact with each other, then the output 
intensity of the aggregate expressed in firings per unit time per neuron 
would be identical with the output of each of its neurons. If, however, 
the parameters are distributed in accordance with certain probability 
frequencies, the output of the aggregate will be an average with respect to 
such distributions. For large aggregates, such an average can be expressed 
as an integral. For example, let the frequency distributions of the refrac- 
tory period 6, the period of latent addition o, and the threshold h (cf. 
Rapoport, 1950b) be given by D(6), S(a), and H(h) respectively. Then, 
if y(«, 6, o, 2) is the output function of a neuron with parameters 4, o, 
and h, and subjected to the input «, and if these parameters are independ- 
ently distributed, we shall have for the output intensity of the aggregate 
(firings per unit time per neuron): 


V (x) =f" ff y@ae, PY DS Sto) H (hy dédodh. (1) 
0 0 0 . 
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If the function y and the frequency distributions D, S, and H are known, 
Y can be computed. It is more likely, however, that experimental pro- 
cedure may allow us to determine Y(x). In that case, although the func- 
tions of the integrand are by no means uniquely determined, nevertheless 
some conclusions can be drawn concerning them on the basis of certain 
assumptions. If, for example, certain forms are assumed for three of the 
functions under the integrals, the fourth may be determined as a solution 
of an integral equation, provided certain conditions which guarantee the 
existence of the solution are satisfied. 

Thus, formally we are dealing with two classes of problems, the direct 
problems, where the output function of the individual neuron and the 
parameter distributions are given, and the output of the aggregate is to 
be deduced; and the inverse problems, where the output of the aggregate 
plus some of the integrand functions are given, and the remaining inte- 
grand functions are to be deduced. 

In the present paper we shall solve some very special cases of the direct 
problem, not so much for their own importance as in an attempt to gain 
insight into some more general aspects. 

Let each neuron of an aggregate be a “‘simple counter,” i.e., it responds 
to every stimulus it receives, provided the stimulus does not fall within 
its refractory period 6. Such neurons were discussed in earlier papers (e.g., 
Rapoport, 1950a). Then the individual output function is 


y (4,6) =xCor4 1)", (2) 
and the aggregate output density will be given by 
oF x D (6) 
Y(e)= fo See. (3) 
For very large frequencies, we shall have 
Viele oF (4) 
0 


We shall suppose that D(6) is continuous to the right throughout, i.e., 
pt ee a fay. (5) 
Since D(5) is a probability frequency distribution, we must have 


J D(a) do=1, 


Hence the only questions concerning the convergence of the integral (4) 
arise in connection with the possible discontinuities of D(5)/5 at the 
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origin. Such questions are not of interest here. We may, however, for 
completeness of discussion, assign the value ~ to Y() if (4) diverges. 
Such may be the case, for example, if our aggregate contains neurons 
with arbitrarily small refractory periods (e.g., if Lim D(O) #0). 


Let us examine the behavior of Y(«) for very large x for two very 
simple distributions of 6, namely, the rectangular and the triangular. 
For a rectangular distribution, we have 


D(6) =0, for OS56<A-<a; 
D (6) = (2¢)—, forA—aSi<A+a; (6) 
D(6) =0,forézA+a. 


Here A is the mean value of 6. 
In this case, equation (3) becomes 


a 4 
Vis)=ae f edb 1, (A+a)x+1 (7) 


Ube es bead 2a A oe i 
and 
toa A+a 


As a approaches zero, 1.e., if the values of 6 are strongly concentrated 
around their mean value, Y() approaches 1/6, as, of course, should be 
the case. As a approaches A, i.e., the rectangular distribution attains its 
maximum spread about the mean value, Y(~) increases without bound, 
since the integral (4) must diverge in that case. Equation (8) may be 
conveniently represented by a series of the form 


De aE a0 Vous 28 ONG 
ei iGs 
1 ia aN 
rs a)teaye |: 
from which it appears that Y() grows monotonically without bound at 
an increasing rate as the spread (a/A) of the distribution increases to 


unity. 
Consider now the triangular distribution given by 


(9) 


Ds) =0, for 6 <A—<¢5 
D(6) =m(é6—A+a),forA—aSb<A; 
.D(8) = —m(i—A-—a), forASi<A+a, 


(10) 


D(6) =0, for 6szA+¢e. 
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The normalization of D(5) requires that m = 1/a’. Therefore 


F 1 4 (§6—A+a)x Ata— (§—A— a)% 3. 
TO aly eae dot f, bx + 1 


a (11) 
1 Ax tax+1 ( >) (Ax + 1)2— a*x? 
=—| ont ae o (Ax + 1)? Ih 
By the application of L’Hospital’s rule, it can be readily seen that 
Limey (a0) ay (12) 
a—0 A 


For the maximum spread, on the other hand, that is, for a = A, we 
have 


2 1 
V Cie ~l4 In 2Aeee +(a+ =)in mai), (13) 
so that 
ine 4 . 
Lim ¥ (x) = ieee (14) 


Thus the asymptotic frequency of an aggregate of simple counter neu- 
rons, whose refractory period is spread in a symmetric triangular distribu- 
tion about the average, is modified by a factor of about 1.4 relatively to 
the asymptotic frequency of its average neuron if the triangular distri- 
bution has the maximum spread. We shall see that for any smaller spread, 
the correction factor is smaller, as is intuitively evident. In fact, we shall 
derive a series, giving the correction factor in terms of a and A. 

Note that setting « = © in equation (11) gives us the asymptotic fre- 
quency of the triangularly distributed aggregate, 


reo) ALLS aaa faa D aa] 
13} 
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[(A+a)lIn(A+ a) + (A—a)In(A—a) —2AlnA]. 
This may also be writtten as 


4 jam(1-2)+efa(42)-m(1-D]}- a 


Since a < A, we may expand the logarithms obtaining 


al ~-a .)+20 (G+ gat.. aE ap 
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and, upon combining terms, 


Y («) =<(1+25+....), (18) 


which is a series of positive terms only and hence monotone increasing 
with a. 

Consider now any unimodal distribution D(5) symmetric about its 
mode (= mean). Let the value of 6 be more concentrated around the 
mean A than those of a triangular distribution with the same range and 
mean. Then we shall say that the distribution D(6) is “sharper” than the 
corresponding triangular distribution. More precisely, if, when we super- 
impose the two distributions one over the other, the ordinates of the tri- 
angular distribution are greater at the extremes and less near the mean 
than those of the distribution in question, then we shall say that the dis- 


DS) 


J 


FIGURE 1 


tribution in question is sharper than the triangular. Such a distribution 
is seen in Figure 1. 

Since for symmetric distributions the area underneath the curve to the 
left of the mean is always equal to $, it follows that the two shaded areas 
in Figure 1 must be equal. But the contribution of the horizontally 
shaded area to the total output is greater than that of the vertically 
shaded area, because the neurons contributing to the former have smaller 
refractory periods. Hence on the left of the mean the total output of the 
triangular distribution is greater. The situation is reversed on the right 
side of the mean. However, the effect of the left side predominates, again 
because of the smaller refractory periods. 

The asymptotic output of an aggregate with any symmetric distribu- 
tion of refractory periods will always be greater than that of an aggregate 
all of whose neurons have a refractory period equal to the mean of the dis- 
tribution. This again follows from the preponderance of the contribution 
to the aggregate output by the neurons whose refractory periods are 
smaller than the mean. Thus we have the result that in any symmetric 
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unimodal distribution sharper than the triangular one of the maximum 
possible spread (a = A), the asymptotic frequency of the aggregate must 
lie between 1/A and 21n 2/A ~ 1.4/A. For more accurate estimates, one 
may compare the distribution with a triangular one of the same range. 
The asymptotic output of the latter is given by the series (18), and the 
asymptotic output of an aggregate with a sharper distribution will not 
exceed it. 

Thus the behavior of the input-output curves of aggregates with sym- 
metric “bell-shaped” distributions (sharper than triangular ones) is not 
appreciably different from that of a single neuron with refractory period 
equal to the mean of the distribution. It can be readily verified that any 
such curve will have a derivative equal to unity at the origin and will 
approach its asymptotic value. Since the asymptotic values will not differ 
markedly from that of the single ‘“‘average’’ neuron, we see that in general 
the input-output curves of most aggregates with quasi-normal distribu- 
tions can be conveniently approximated by the input-output curve of the 
average representative neuron. 

The Smoothing-oui Effect of Distributions. If the stimuli received by the 
neurons in an aggregate are regularly spaced in time instead of Poisson- 
distributed, a somewhat different problem arises. It was shown earlier 
(Rapoport, 1950a) that the input-output curve of a single neuron sub- 
jected to regularly spaced stimuli is discontinuous, namely (cf. Rapoport, 
1950a, Fig. 1), 


y (#,6) = %, for OSx% el. 


Oe 
SCRE) S for Su <5; 
6 6 (19) 
x — k 
6) == ——s =. 
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The output of an aggregate of such neurons, whose refractory periods 
have a distribution D(6) will then be 


V(x) = fy (#8) D(a) do. (20) 


We can now ask what will be the “smoothing-out” effect of the dis- 
tribution function D(6), in particular, what conditions must be satisfied by 
D(6) in order that V(x) be monotone increasing with «? 

As in our earlier paper (Rapoport, 1950a), we write (x, 6) as a func- 
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tion of 6, holding x fixed, namely, 


W126) = 0 for 0S 8 <=; 


y («, 6) =5 fore se 
: “ (21) 
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Then equation (20) becomes 


/z 2/x 
Y (x) =a(f- D(o)do+5 f D(s) ds 
etd 1 kz 
+....g f D(ads). 


If D(6) is continuous in 6, Y(x) is a continuous function of x. We seek 
the condition on D(6) such that Y(x) is a monotone increasing function 
of x. We shall assume that Y(x) possesses a derivative everywhere and 
ask under what conditions Y’(x) > 0 for all x. 

Denoting 1/x by z, we re-write equation (22) as follows: 


Z(2) =2(f D(a da+3 f D(d)dd+....). (23) 


The condition Y’(x) = 0 will now be replaced by the condition 
Z'(z) < 0. Differentiating (23) with respect to z, we obtain 


a (a) =-4 (fp (a) dst4 f D(a) dd+....)+F1D(®) on 


=—D(0)+D(22) —4D(z) + D (32) —3D (22) +....). 


We shall assume that Y(~) < o and hence that D(0) = 0. Then, col- 
lecting terms in (24), we obtain 


Boi ia —A(f diaast+3f D(a ast. ps .) 


zs 
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(25) 


4+=[}D (2) +4D (22) +3D (32) +....1. 


The expressions in both brackets are positive, so that to have Z’(z) < 0 
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for all z, we must have 


z 22 
Male Gog eee 
(foo aasta f (8) ds+ ) ie 
7D Ge. 


for ail z. Inequality (26) is the condition on D(6) sought. 

We shall investigate some special cases of D(6) with regard to their 
“smoothing-out” effect on the aggregate output under the condition of 
regularly spaced stimulation. Consider, in particular, the rectangular 
distribution given by equations (6). Applying equation (22), we see that 
so long as « is sufficiently small (i.e., 1/a > A + a), V(x) = x. Now let x 
increase and let us see what happens as 1/x just passes (A + a) from right 
to left, i.e., becomes less than that largest value of 6 in the aggregate. 

We note that d(1/x) = —d«/x?. In particular, when x = (A+ a)" 
d(1/x) = -—dx(A + a). Since 


At+a 
{iD C3) ds ae 
0 


the expression within the parentheses of the right side of (22) becomes 


(fd a do43 ff “D(a dd). (27) 


(Note: We have broken off the series within the parentheses at the second 
term, because we are computing V(x) for those values of x for which 
1/e <A+a < 2/x.) Let x = (A+a)++ dx. The expression within 
the parentheses of (27) then becomes 


,_ ax Ara? 


a (28) 


where the quantity subtracted from unity represents the loss of fre- 
quency due to the fact that the neurons for which the input frequency 
has just passed the critical frequency (A + a)— are responding to only 
every other stimulus. The total frequency per neuron for « = (A + a) + 


dx is thus 


[ (A+) 4dx](1- SEEM") (29) 


Neglecting higher order infinitesimals, we may write it 


(A+a) 4d (1-22), (30) 


Hence, at the critical input frequency (A + a)- the output curve will 
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remain increasing if, and only if, 
A+a< 4a; ee (31) 
Inequality (31) gives the minimum spread of a rectangular distribution 
which will “‘smooth-out” (in the sense of introducing monotonicity) the 
output curve at the first discontinuity. The situation at the succeeding 
discontinuities can be similarly computed. 
The derivative Y’(«) = 1 in the range 0 S x < (A+ a)—. The com- 
putation of V(x) at the points ae ae to the right of (A + a) gives 


3a — 
Va ae =a (32) 
Therefore Y’(«) in that region is given by 
: 3a — 
Yes) ——. (33) 


We thus see that there will always be an abrupt change of slope in the 
output curve at the first critical point. The rectangular distribution can- 
not smooth out the output curve to make the derivative continuous. We 
shall see, however, that the triangular distribution can produce such a 
smoothing-out effect. For this case, if A+ a > 1/x > A, equation (22) 
becomes 


1 1 1/z 1 At+a ) 
cea wee y ($—A—a) dd—7 f GA ards 


2 @ 
(34) 
AG (oe en ee 
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Setting x = (A+ a)-", we obtain 
VilA+e)J=1. (36) 


Thus not only is there no decrease in output initiated at the first critical 
frequency, but the output curve is continuous over it with a continuous 
derivative. 

Nevertheless, Y(«) may or may not be monotone if D(6) is a triangular 
distribution. We shall make a sample investigation in the region where 
(A + a) <« < A“ to determine under what conditions Y’(x) can 
change sign in that interval, i.e., under what conditions Y(«) will fail to 
be monotone. For this to happen, we evidently must have holding 
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simultaneously 


3 1 A? A 
= = =3{))< Sul 
Ne Ee 4a? 2a ’ om 
(Aid) tS ec Noe (38) 


Moreover, since D(0) = 0, it follows that A > a. Hence A+ a < 24. 
This condition insures that equation (27) describes Y(*) in the entire 
range of « given by (38), since as 1/x approaches A from above, 2/x has 
not yet passed A + a. 

Solving (37) for «?, we obtain 


42 = (A?+ 2aA— 3a?) 7. (39) 


Substituting this value into (38), we see that the following inequalities 
must be simultaneously satisfied: 

A?+ 2aA— 302s (A+a)?; (40) 

ME 2ah — Sat (41) 

Inequality (40) is identically satisfied. Inequality (41) holds if, and 

only if, 
A 
a<2z. (42) 


Thus inequality (42) is the condition that Y («) fails to be monotone in the 
interval (A + a)-! < « < At, It follows that if Y(x) is to be monotone 
in that interval, we must have 


CHP. (43) 


wl e& 


which gives the minimum spread of the triangular distribution insuring 
monotonicity in the region (A+ a) S x < At. 

Proceeding in this way, we can examine the behavior of Y(x) in suc- 
ceeding intervals to see if more drastic conditions are required for its be- 
ing monotone. One suspects that (31) is sufficient for the entire range in 
the case of the rectangular distribution, and, similarly, (43) is sufficient 
for the triangular distribution, since the satisfaction of these inequalities 
removes the greaiest oscillation of the discontinuous function (x, 5), 
namely, the first one (at x = 1/6). The proof of this conjecture is per- 
haps not difficult. 

The general qualitative result one may deduce from the above discus- 
sion is that given a sufficient dispersal of 6 (measured by a in the triangu- 
lar distribution), the output curve of an aggregate of neurons subjected 
to regularly spaced stimuli can be smoothed out. Similar results are likely 
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to be obtained from more general distributions. In general the smoothing- 
out effect will depend on the dispersal of the distribution. However, the 
converse of this conclusion is also interesting. If the dispersal of a distribu- 
tion is small, the aggregate may exhibit a smooth output curve in re- 
sponse to regularly spaced stimuli, and yet the curve may not be monotone 
(i.e., exhibit maxima). Here then is a theoretical possibility of having non- 
monotone output curves without postulating any inhibitory action 
whatsoever. 

Now a regularly spaced sequence of stimuli can be considered as a dis- 
tribution of intervals without dispersal. In general a time sequence of 
stimuli will not be absolutely regular, so that the distribution of the in- 
tervals between the stimuli will exhibit a finite dispersal. We have seen 
that the dispersal associated with a Poisson shower is sufficient to give a 
monotone output curve even of a single neuron (or an aggregate with zero 
dispersal of refractory periods). On the other hand, if the dispersal of the 
intervals between stimuli is zero, certain minimum dispersal of refractory 
periods is required for monotonicity. It follows that monotonicity can be 
obtained with some sort of combination of interval dispersal and refrac- 
tory period dispersal. If this combined dispersal falls below a certain value, 
non-monotonicity appears. We have seen, however (Rapoport, 1950a), 
how a non-monotone output curve can be used to obtain sensitivity to 
particular ranges of frequency. It now appears that such an effect can be 
obtained by a “tightening up” of either the firing pattern or of the dis- 
tribution of refractory periods or both. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Labora- 
tories and The University of Chicago. 
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In the theory of diffusion of electrolytes the following assumptions are fre- 
quently made: (i) the electrolytic solution is electrically neutral everywhere, (ii) 
the ionic concentrations and the electric potential all depend on a single Carte- 
sian coordinate as the only space variable. Often the electric potential of the solu- 
tion is determined on the basis of the Poisson equation alone, disregarding any 
other relation between this potential and the ionic concentrations. Since the 
Poisson equation only represents a condition which the potential fulfills, the use 
of this equation alone may lead to error unless the explicit relation for the po- 
tential involving a space integration of ionic concentrations is also taken into 
account. But if this relation is used the Poisson equation becomes redundant 
and, more important, assumptions (i) and (ii) appear unacceptable, the former 
because it leads to a zero electric potential everywhere, the latter because it is 
mathematically incorrect. The present paper is based on general equations of dif- 
fusion of ions, excluding the Poisson equation. These equations form a system of 
nonlinear integrodifferential equations whose number equals the number of ionic 
species present in the solution. It appears that when allions are distributed sym- 
metrically around a point all functions related to the above system of equations 
can be made dependent on a single space coordinate: the distance from the center 
of symmetry. Two methods of successive approximations are given for the solu- 
tion of the equations in the case of spherical symmetry with limitation to the 
steady state. These methods are then applied to the study of the distribution of 
ionic concentrations and electrical potentials inside a cell of spherical shape in 
equilibrium with its surroundings. These methods are rapidly convergent; exact 
theoretical values of the electric potential are calculable on the boundary of the 
cell. It appears that the potential at the center of the cell is not more than 
~50% higher than at its boundary and that variation of concentration in- 
side the cell is not very large. For instance, with 100 mV on the boundary the 
ionic concentration there is about four times higher than at the center. Calcula- 
tions show that extremely small amounts of electricity are sufficient to account 
for the electric potentials currently observed. In a cell of 100 micra diameter an 
average concentration of only 10-4 mole/cm# of a monovalent ion would be suffi- 
cient to give 1 millivolt on the boundary. This concentration is directly propor- 
tional to the voltage and inversely proportional to the square of the cell diameter. 
Most of the numerical results given above are obtained by considering only 
those ions whose electrical charge is not compensated for by ions of an opposite 
sign. The total concentrations may be much higher than those quoted. The 
theory does not take into account possible effects of structural heterogeneities 
which may exist in the cell, particularly of various phase boundaries. An inciden- 
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tal result shows that the Boltzmann distribution function in the form employed 
in modern theory of electrolytes is fundamentally a consequence of the mathe- 
matical theory of diffusion alone. It is pointed out, however, that Boltzmann dis- 
tribution is not always compatible with the definition of the electric potential. 


I. Some general equations of diffusion of ions. The classical theory of 
diffusion of electrolytes is based on the differential equation 


6) Ci ae - ; ae c | 1 
er D; | divgrad Citar div (c; grad f) |. (1) 
This equation is an expression of the fact that ions move under the action 
of a concentration gradient grad c; and of the electric potential f, which is 
the sum of an externally applied potential, if any, and of the potential 


1 ¢:(Q) dVoa 
f(P) =E LS, RE, ae (2) 


due to the ions themselves. In the above expressions ¢;(Q) is the concen- 
tration of the ions of kind z at Q; an infinitesimal volume element around 
this point being indicated by dVo; f(P) is the potential at P, R(P, Q) is 
the distance between P and Q, the charge of the ion of kind 7 is indicated 
by ¢,, and the dielectric constant of the solution by K, the diffusion coeff- 
cient of the 7th ion in the solvent is indicated by Dj, the interaction be- 
tween the ions of various kinds is assumed to be of a purely electrical na- 
ture, that is, mutual diffusion effects are neglected. Electrical forces be- 
tween the ions of the solute and the polarized molecules of the solvent are 
also neglected. The meaning of the other symbols is: ¢ = time, k = Boltz- 
mann constant, and T = absolute temperature. The mobility of the ions 
is taken to be u; = D;/kT, which implies an application of the law of 
perfect gases to the ions in solution. 

We will limit ourselves in this note to the case in which no external elec- 
tric field is acting on the solution. Under this condition equation (1) was 
used by M. Planck (1850) in his theory of electrolytic solutions. Instead 
of equation (2) he used the Poisson equation 


€; Cy 


divgrad f= — 40 >), (3) 


but neglected its right-hand side. K. Sitte (1934, 1935) pointed out that 
this procedure is not justified because, using equation (2), it would be 
equivalent to putting electric potential equal to zero everywhere, which 
would reduce equation (1) to a diffusion equation of neutral molecules 
and not of ions. J. J. Hermans (1936) tried unsuccessfully to justify this 
approximation by a mathematical argument. Of course, a method could 


DIFFUSION OF ELECTROLYTES 87 


be devised for the solution of the above equations by expanding the func- 
tions in MacLaurin series of 1/K. Such a method has been used in extend- 
ing the Debye-Hiickel theory of electrolytic solutions (Gronwall, La Mer 
and Sandved, 1928; La Mer, Gronwall and Greiff, 1931). As in Planck’s 
theory the zero order terms in such a method would imply a zero right- 
hand side in equation (3). 

If the concentrations c,(i = 1, 2,...,) are considered as the un- 
knowns, the functions c; are solutions of 7 nonlinear integrodifferential 
equations obtainable by substituting f from equation (2) into (1). Instead 
of these equations Planck and Sitte worked with equations (1) and (3) 
and made the additional assumption that all the functions c; and f depend 
on one Cartesian coordinate as the only space variable. However an in- 
spection of equation (2) shows that the electric potential f and the density 
of the electric charges €;c; to which f is due cannot all be functions of only 
a single Cartesian coordinate. 

Since equations (1) and (2) are very difficult it is worthwhile to look for 
cases in which concentrations and electric potential can be made dependent 
on a single space coordinate. We will show that this is the case if the 
system has spherical symmetry. If the ionic concentrations c; depend 
only on the distance r from a fixed point we have 


f= Ze foredrt f° redr), (4) 


where r = 7; and r = 72 are two bounding surfaces of the solution with 
ry <1. If r, = 0 the solution fills a sphere of radius 72. Writing equa- 
tion (1) in spherical coordinates we obtain: 
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and eliminating f we have 
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where k; = 42e€;/KkT and 


Q=4r ef reds (7) 


a) 


is the total electric charge existing within a sphere of radius 7 and center 
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at the origin. Equation (6) can also be written as: 


r? dc;_ 9 ; 2 Set) he eee. (8) 
D; Ot or 


where h; = ¢,/KkT. It is seen from equation (7) that if the electric charge 
Q(r) in a sphere r = const. is a maximum or a minimum the surface 7 = 
const. is electrically neutral. From equations (4) and (7) we have for the 


electrical potential: 
QO An "2 ; 
a + : af rco;dr. (9) 


The above equations describe a very general type of diffusion phenomena 
in which c; and f both are expressible in terms of a single space coordi- 
nate r. 

II. The steady state. If sources and sinks of ions have appropriate asymp- 
totic characteristics as # > ~ or if they cease to act ast — o, a steady 
state is reached: lim 0c;/0t—>0 when t— ~, i.e., using equation (5) 
we have: a , 

= €; Jf 

Stes a As (10) 

where the A; are constants. The number of 7th ions flowing per unit time 
at t= © through a spherical surface r = const. is 47D;A;. Equation 


(10) can be rewritten, using equation (8), in the form: 


» dCs 
dr 


eee (11) 


The steady state is ruled by x of these relations, if 2 is the number of ionic 
species. Eliminating Q between any two of these relations we see that 
any c; can be expressed in terms of any other one: 


Bi (rs 16561) = 0564) +a f(r YZ, (12) 


where r* is a value of r at which all the concentrations are known a priori, 
and 


Ex(r*,r) = 7 exp (af, lam). 


Thus the determination of the relation between the various ionic con- 
centrations in steady state requires, in general, not only the knowledge 
of the amounts of ions flowing in steady state through a boundary (con- 
stants A;), but also the knowledge of all the steady state concentrations 
[c:(r*)] at least at one point (r = r*) of the solution. In general these 
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data can be obtained through an analysis of the non-steady state, which 
we plan to discuss in a successive note. There are, however, important cases 
in which these constants can be secured without solving the non-steady 
state equations. We will consider one of these cases in the following sec- 
tion. In any event the system of equations (11) can be changed by 
means of equations (7) and (12) into 2 mutually independent equations 
in the c;’s. These integrodifferential equations are very complicated how- 
ever. 

It can be seen from equation (11) that a uniform distribution of ions in 
steady state is impossible unless the solution is everywhere neutral, that is, 
the electric potential is zero everywhere. In fact the condition that c; be inde- 
pendent of r implies Q = constant, which is possible only if Q is identically 
zero. 

From now on we will limit ourselves to the case in which there are 
no sources and sinks of tons acting att = ~. Then A; = 0 and equation 
(10) gives 


Gs—> €3(0) €s CF) where e;(f) =exp ae (13) 


and ¢,9) are constants depending only on 7. The above is the well known 
Boltzmann equation in the form currently used in the Debye-Hiickel theory 
of electrolytes (see e.g. Harned and Owen, 1950, p. 25). It is interesting 
that this equation can be derived from the mathematical theory of diffusion 
alone which is based on no assumptions beyond the principle of mass con- 
servation, Fick’s law of diffusion, the formula u; = D,/kT implying the 
law of perfect gases, and the assumption of direct proportionality of ionic 
velocities to electric forces. However the Boltzmann equation is not always 
compatible with the definition of the electric potential. In fact if the concen- 
trations c; are all functions of one Cartesian coordinate only, equation (13) 
requires that the potential f also be such a function. But this contradicts 
equation (2). Thus, to a thermodynamical inconsistency of the Boltzmann 
equation noted by L. Onsager (1933), a new one is herewith added. 

If there is a point in the solution at which the potential is zero, the 
concentrations c; at that point equal c,o). If N; is the total number of 
ions of kind 7 existing at ¢ = ~, then 


4nci) f e:(f) dr = Ni. (14) 


From equation (13) we have the following relation between concentra- 
tions and electric potentials at any two points r = a andr = b: 


fa) —f (0) =n (15) 


90 I. OPATOWSKI 


The determination of f in terms of 7 requires the solution of the integral 


equation in f: 
“e 1 r ’ tis : 
EI = Decwl(s frenart f re(fdr) (16) 


which is obtainable from equations (4) and (13). 

Equation (16) determines f as a function of 7 and of the ¢,o’s. If N; in- 
stead of cio) are known a priori, we can use the following equation in- 
stead of (16): 


f(r) =o bait Cf reat foreNar), (17) 


i re; (f) dr 


i 


which is easily obtainable from equations (14) and (16). Both equations 
(16) and (17) are convenient for a determination of f(r) by successive ap- 
proximations. The potential on the external boundary of the solution 
is easily calculable as follows: 

SD: e,N 4 


f(r) = Poy (18) 
This formula, which is also a simple direct consequence of the definition 
of the electrostatic potential, shows that the potential on the external 
boundary (7 = re) is zero if the solution is totally neutral 


(Sani=0), 
For the potential of the internal boundary we have 


Sani fre (Nar 


a 


f(n) = a3 : (19) 
K fre (far 
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III. Solution of the equations of the steady state. Electric potential and 
tonic concentration in a spherical cell. 

A. Method 1. Consider as before a solution bounded by two spherical 
surfaces of radii y = 7, and r = fz, or a solution filling out a whole sphere 
of radius rz. If f does not vary a great deal in the solution we can put f = 
constant as a starting approximation in an iteration process based on 
equation (17). We then obtain an approximate expression fi(r) of the 
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potential f(r) as follows: 


3 


1/3n-72 , i 
Ai) =z 7-4). (20) 


ys v3 
if 's nm 


Since |/f,(7)| is a decreasing function of 7, thus—in a first approximation 
at least—the electric potential has its maximum magnitude at the center and 
minimum on the boundary. At r = 7, we have 


(rit re) Des 


A= (21) 


oF (Peee eee et 
Therefore, the steady state potential at the center of the sphere (r; = 0) 
equals, in first approximation, three halves of its value at the boundary. By 
putting f = f, in the right-hand side of equation (17) a second approxi- 
mation, fo of f, is obtained. It can be expressed in terms of tabulated 
functions if the solution fills a whole sphere, that is, if 7: = 0, 


‘ exp (He; rs) + exp (He; 71?) dr 
ro 
fo (7) aK. eae ae ee) 
aoe exp (He; 72) -—f exp (He; r?) dr 
ro 0 


where 

ae e;N Z 

sit iMG) 
TDK Der 

[For a bibliography of the tables of the integral appearing in formula (22) 
see Fletcher, Miller and Rosenhead, 1946, p. 219.] Both approximations 
fi and fz give the exact expression of the potential f atr = rz. At r = 0 the 
approximation f2 gives 


fa (0) = ———— exp Heir) — Nels (23) 
ERE exp (He; rs) — (r2V He) malf “exp (x?) dx 
0 


As an application let us consider a spherical cell of 100u diameter in 
steady state, with no ionic sources inside, and in equilibrium with its sur- 
roundings. A part of the positive and negative ions existing in the cell 
neutralize themselves. Therefore they can be disregarded in studying the 
electrical properties of the cell; these are due only to the ions giving a net 
contribution to the density of the electric charge in the cell. For simplicity 
we assume that these latter ions are all of the same valence and call e the 
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charge of each of them and N their total number in the cell. The net elec- 
tric charge of the cell is then We. Since the sign of the charge does not 
affect the formulae except for their sign, let us assume Ne > 0 and the 
corresponding ions as monovalent. If we put K = 80 then formula (18) 
gives us the number of ions in question per each mV of electric potential 
on the boundary of the cell, as follows: 

N Kr 


-= — = 8 i . 
Fra) = 300,000. 2,780 ions/mV 


Thus less than 3,000 ions are sufficient to produce 1 mV voltage on the boundary 
of a cell of 100u diameter. This number of ions is directly proportional to 
the voltage and to the cell diameter. Since the dielectric constant of a cell in- 
terior is likely to be smaller than 80, the number of ions required is also 
likely to be smaller than the above figure. In any event this is a very small 
number of ions. This fact appears particularly clear if we consider the 
average ionic concentration in the cell Cay which gives 1 mV on tts boundary 
(monovalent ions): 


CaS 2 ioe 1o-u mole/cm? mV . 
2 


With K = 80and 27, = 100p this gives ca, = 8.8 X 10-% mole/cm’ mV, 
a quantity inversely proportional to the square of the cell diameter. 

Both Hr} = f(r2)/2kT and f2(0)/f,(0) depend on the potential f(r2) at 
the cell surface but not on its radius 72. If f2(0) ~ fi(0), fir) is a good 
approximation to f(7). This fact is seen to depend on f(r2) but not on 7». 

The table below gives the first and the second approximations to the volt- 
age at the center of the cell {,(0) and f.(0) for various voltages f(r) on the 
boundary of the cell at T = 300° K. The results do not depend on the cell 
diameter. 

f (re) in mV 1 50 100 


fi(0)in mV ERS i 150 
fo(0)in mV ie 69 128 


It is seen that the approximation method for the calculation of f is quite 
good, particularly for lower voltages. For higher voltages the result can 
be improved by calculating a third approximation. This is obtained by 
putting f = f: in the right-hand side of equation (17) or (19) and carry- 
ing out numerically the necessary integrations. As the voltage f(r2) on the 
external surface of the cell tends to zero, Her? +0 and equation (23) 


gives f2(0) — 3f(r2)/2, in agreement with the first column of the above 
table. 
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B. Method 2. In the previous successive approximation method the 
starting approximation was a constant potential throughout the cell. In- 
stead of the latter we can take a suitable step-function as the starting ap- 
proximation. We illustrate this procedure by calculating the potential of 
the cell using the data of the previous example. As a starting approxima- 
tion we put f(r) = f(0) for0 <r < 72/2 and f(r) = f(re) for r2/2 <7 < rr. 
Since f(72) is directly calculable from equation (18), relation (19) gives a 


transcendental equation for f(7:): 
3X. — 44 
exp (%1 + %2) = Waeet0 Uc 

where «; = ef(r,;)/RT, with 7; = 0. The solution of this equation for r2 = 
100u, 7: = 0, T = 300° K, f(r2) = 100 mV gives a first approximation to 
the potential at the center f:(0) ~ 135 mV, which is in quite good agree- 
ment with the results of the previous method. Equation (17) can now be 
used to determine fi(7), an approximation to the potential f(r) in the 
whole range 0 <r < ro: 


_ €N a,+ 4 exp (41 — %2) 


Ma RG 1+7 exp(*,—%) ’ 2 (4) 
where 
(3 ra— Ar for rst 
sa rs 2 
aa Salis 
= for "25 
Or; for SZ 
i aes ro 
127,-—-—— 4r? for zope 
7 


If, following method 1, we take x; = a, formula (24) reduces to (20). A 
second approximation to f(r) can be obtained by putting f(r) = f(y) in 
the right-hand side of equation (17) and calculating the corresponding in- 
tegrals. The ratio of local concentrations on the boundary and at the center 
is c(r2)/c(0) ~ 3.9, according to formula (13), in the case of 100 mV on 
boundary, if we assume 135 mV as the correct figure for the potential at 
the center. 

The following modification of method 2 has been suggested to the writer 
by Professor H. D. Landahl. Since, according to equation (17) 


oa eV; 


OP) ee (32 Sate 
or peepee *"\ dar r=r, Tor; : 
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we can use as the starting approximation to f(r) an appropriate function 
having the calculated values f(7:) and f(r2) as well as the above slopes 
df/dr at r = r, and r = 72. Such a function can be easily represented by 
a curve since f(r) has no maxima or minima in 7 <r < 72 [cf. equation 
(17)]. This modification is likely to yield a more rapidly convergent re- 
sult than method 2. It requires numerical integrations throughout, how- 
ever, unless the curve can be fitted by a suitable analytical expression. 

The distribution of potential in a spherical cell has been treated recently 
by F. Booth (1951) using the method mentioned above of Gronwall, 
La Mer and Sandved. Numerical computations of the potential using this 
method are laborious and the power series involved presents several dis- 
advantages which have been discussed by Booth (loc. cit.). 

This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of The University of Chicago. 
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The paper outlines a possible further development of the suggestion made by 
the author in his recent book. Ideology is defined here as a verbalized, or at least 
verbalizable, behavior pattern which may be adopted by society. After a brief dis- 
cussion of possible classifications of ideologies, a study is made of those ideologies 
which refer to the question of social and ethical interrelations. Different kinds of 
ideologies may be represented by behavior matrices, introduced in the author’s 
book. A uniparametric representation of all such matrices is suggested and dis- 
cussed. Next the previous results on social imitation are extended to the case of 
n different behaviors, each of which is determined by a particular value of a con- 
tinuously varying parameter. It is shown that, depending on some other social 
parameters, changes from one ideology to another may proceed either quasi- 
continuously or definitely discontinuously. The paper concludes with some general 
speculations on the possibility of applying the above results to a mathematical 
interpretation of history. 


In our recent book (Rashevsky, 1951a) we suggested the possibility of 
developing a quantitative dynamics of ideas, considering ideas as complex 
verbal behavior patterns. From a neurobiological point of view the overt 
behavior and the covert idea of the behavior merely form two parts of 
the same neurobiophysical process. The word “idea” was understood 
as meaning merely a verbalized, or at least verbalizable, statement of the 
process of behavior. It has been shown by E. Jacobson (1931) that when 
an individual thinks of reciting a given poem action currents are produced 
in the muscles of his tongue which show a time pattern identical with 
that of the action current observed when the poem is actually recited. 
However, the intensity of the action current in the latter case is about 10° 
times greater. This and other experiments by Jacobson (1930) make it 
plausible to assume that the thought or idea of a behavior is the result of 
a neurobiological process which forms an essential part of the neurobio- 
logical processes that result in the behavior itself. Thus no metaphysical 
meaning whatsoever is attached here to the word “‘idea.” 

However, since this word has to many persons a definite connotation, 
different from the one used here (e.g., Platonic ideas), we shall use instead 
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the word “ideology,” which, though perhaps not as appropriate, is free 
from the connotation mentioned above. However, we shall also use the 
word “‘idea’”’ as defined above. 

It is interesting to note that as far back as 1909 Henry Adams (1949) 
suggested the possibility of a dynamics of thoughts or ideologies which 
might have thrown light on the understanding of history. In his sugges- 
tion he attempted to construct a dynamics of ideologies along the pat- 
terns of dynamics of physical systems. His approach somewhat resembles 
Haret’s approach to social dynamics. We have elsewhere (Rashevsky, 
1948, 1951a) pointed out the fundamental weakness of a mere mimicking 
of physical sciences in dealing with social problems. 

In this paper we shall give a very sketchy outline of a dynamics of 
ideologies based on the mathematical biology of social behavior, as de- 
veloped previously (Rashevsky, 1951a). 

I. Quantitative Classification of Ideologies. We shall consider two main 
qualitative subdivisions of ideologies and within each we shall attempt a 
quantitative classification. 

The first subdivision contains the ideologies concerned with the nature 
and causes of things. They may be called causational ideologies and spring 
from attempts to answer the question ‘‘why?” All beliefs in deities of 
various kinds and ideas about their parts in creating and operating the 
universe belong to this category, and form the first subcategory. To the 
same category, but to a different subcategory, also belong different scien- 
tific ideas concerning the nature of the world. The fundamental difference 
between the two subcategories is that in the former the acceptance on 
faith plays a cardinal role, whereas in the latter rational critical analysis 
is paramount. As we have pointed out elsewhere (Rashevsky, 1951a, chap. 
xxi) complete freedom from acceptance on faith is a physical impossibility. 
We cannot verify every statement which we hear and which we intend to 
use. However, whereas some individuals will accept some statements on 
faith as a matter of expediency, and will always admit that the statement 
so accepted may turn out to be wrong, others consider the slightest doubt 
of certain statements as inadmissible. Herein lies the fundamental dif- 
ference between the scientific and religious, or rational and arational (loc. 
cit., p. 170), attitudes. A neurobiological mechanism which explains how 
the same individual may exhibit both attitudes has been discussed (Joc. cit., 
chap. xxi). 

We may call the ideologies of the first subcategory arational, that of the 
second rational. In this paper we shall restrict ourselves to the former. A 
theory of the latter will be developed subsequently. 
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A quantitative classification of that subcategory may be attempted by 
considering the number of different deities involved in the various beliefs. 
We thus find a quasi-continuous change going from the innumerable 
“gods” of the fetishist through the “higher” polytheistic religions to the 
so-called monotheism of Christianity and Islam. Together with the reduc- 
tion of the number of deities goes an increase in their abstraction. The 
process may be readily understood in terms of man’s developing knowl- 
edge of natural phenomena. In the earliest stages of development man 
had not learned the fallacy of false generalization. An accidental coinci- 
dence between two phenomena led him to consider one phenomenon as 
the cause of the other. The two phenomena became permanently asso- 
ciated in his mind in a causal relation. A “‘deification”’ of different objects 
or phenomena thus occurs. With progressing knowledge those original 
associations, which are formed as conditioned reflexes, are gradually dis- 
solved. An individual who learns that he can make a useful tool of which 
he has full command from a piece of wood which he originally deified will 
deify it no longer. Thus the total number D of deities is being reduced 
through knowledge. As man gains more and more control over nature, the 
concept of deity becomes more and more abstract. Man will not deify an 
object or phenomenon over which he has control. 

Thus one of the first relations which may be looked for in a dynamics of 
ideas is a relation between D and the total amount K of knowledge. To 
‘this end we must first find a proper measure of knowledge, and then es- 
tablish the law of its variation with respect to time. Possible expressions 
for K have been tentatively suggested in terms of information theory 
(Rashevsky, 1950). A development of a theory on that basis may be at- 
tempted. 

The amount of knowledge is in general distributed non-uniformly in a 
society. A few individuals have much knowledge, many have little. The 
general theory of social distributions, developed in chapter viii of Joc. cit. 
may be applied to this case. The problem is somewhat simplified by the 
fact that when two individuals who possess different amounts of knowl- 
edge interact, the individual with less knowledge may gain some, while the 
individual with more knowledge cannot lose any. Thus the function 
p(s, s’, A), where s and s’ now denote the amount of knowledge, is such 
that for s-> s’, A= 0, while for's < s,A > 0. 

Due to a process of imitation, ‘such as discussed in chapter xiii of Joc. 
cit., it is possible for a social group to adopt, so to speak “officially,” an 
ideology which corresponds to the higher knowledge of the few, although 
the majority may have a much lesser knowledge. The less enlightened ma- 


98 N. RASHEVSKY 


jority, though adopting the ideology of the more enlightened minority, will 
likely exhibit some features of behavior which correspond to the lesser 
level of knowledge. Thus in some Christian countries beliefs in ghosts, 
mermaids, etc. still persevere among non-educated people. 

The opposite effect, however, may also occur: the society as a whole 
may adopt an idea which corresponds to the lesser knowledge. 

As in other cases of imitative behavior, changes from one set of ideolo- 
gies to another will in general occur discontinuously, even though K may 
vary continuously. 

In a preceding paper (Rashevsky, 1951e, hereafter referred to as Beliefs) 
we suggested a neurobiophysical mechanism for the formation of various 
beliefs. According to that picture the adoption or non-adoption of some 
religious belief by a society is determined by the interaction of the process 
of early conditioning to the belief with the growing tendency to ana- 
lyze everything critically. The strength of conditioning caused by ex- 
ternal factors is determined by the quantity b in equation (2) of Beliefs. 

Now consider a society composed of 7 spatially distinct subsocieties 
which closely interact with each other. Let the belief patterns to which the 
individuals are conditioned in general be different for the different sub- 
societies. However, let all those patterns contain one common element. 
Thus the pattern for the kth subsociety may be represented as consisting 
of the set of elements 


(1) (2) 


Bee BeEB Rew Be: (1) 


where B{" is in general different from Bie} and where Bo is the common 
element. Due to the interaction of the subsocieties an individual in the kth 
subsociety will be principally conditioned to the pattern B;, but to a lesser 
extent he will be conditioned to the remaining (m — 1) patterns B;(i ¥ k). 
Let the relative frequency of conditioning to any other pattern be ex- 
pressed by the fraction w(< 1), if the frequency of conditioning to B, is 
taken as unity. The greater the interaction of the subsocieties, the 
greater w. If w = 1, this means such a close interaction that each individu- 
al is equally exposed to the conditioning to the beliefs of his own sub- 
society and of any other. 

The common element By is conditioned every time any pattern is con- 
ditioned. Other elements are conditioned only when their corresponding 
pattern is being conditioned. Assuming approximate linearity between 
the final strength of conditioning and the number of repetitions, we find 
that the strength of conditioning of By will be 1 + w(m — 1) times greater 
than the strength of conditioning of any other element B{"). 
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In line with the developments presented in Beliefs, if we now ask what 
the probability is that any of the partial beliefs B{” will be abandoned by 
the society as a whole, we find that this probability depends on the values 
of the constant b (cf. Beliefs) for each of these beliefs. From what we have 
just said it follows that the value 6 of for the belief Bo is 1 + w(n — 1) 
times as large as the value d§” of b for any of the other beliefs B{”. Thus 
we have 


b = Tlf o(n—1)1 be. (2) 


In absence of imitation the fraction pv of individuals who do not accept Bf” 
is given by equation (49) of Beliefs, into which equation (20) should be 
introduced for ¢,,. This gives 


bey Mae Xr een 
| (2 psa Z| 
eS Bags E 
i 2 Pt: av s?+ 6% i 3) 


If \/¥ is sufficiently large, then for sufficiently small values of b the 
fraction v will approach unity. For sufficiently large values of w and n, 
the values of b{") may be sufficiently small to insure a non-acceptance of 
the beliefs; yet the value by may be, because of (2), still so high that be- 
lief Bo is accepted by the society as a whole. 

Thus the condition for the survival of the common element of several 
belief patterns is a sufficiently large number x of subsocieties and a suff:- 
ciently strong cultural interaction w between these subsocieties. 

Such conditions are realized by large empires composed of numerous 
countries with different belief patterns. Let B; stand for any polytheistic 
religion. The common element By is the belief in the existence of deity in 
general. The number of deities and their attributes constitute the ele- 
ments BS”. 

The above discussion may throw some light upon the question of why 
Christianity, which is relatively monotheistic, spread approximately at 
the time of the largest extension of the Roman Empire. Other ancient em- 
pires were also composed of several subsocieties, but the value of m, and 
possibly of w, was the largest for the Roman Empire. 

We may also ask the question: What are the characteristics necessary 
to the particular subsociety in which the monotheistic idea first origi- 
nated? Clearly, that subsociety will have values of b{? which are already 
rather small compared to b. Such characteristics will be present in a 
subsociety which previously has had cultural contacts with numerous 
other societies which have different belief patterns. This condition may 


100 N. RASHEVSKY 


be fulfilled by some nomadic societies, or societies that frequently in the 
past have been temporarily subjugated by others. 

Thus it may not have been a mere accident that the earliest concepts 
of monotheism appeared in the originally nomadic Jewish society, which 
formed a part of the Roman Empire at the time when the conditions for 
a wide spread of monotheism became ripe. 

We now pass to the second important category of ideologies. It con- 
tains ideologies which deal with the question: How shall a society or in- 
dividual behave? This is essentially an ethical question. There seems to 
have been a definite correlation between the ethical codes and religious 
views in the past. This is probably due to the fact that both depend on 
some common variables—perhaps, for example, the knowledge K. 

Quite generally the types of behavior of a society may be described by 
the behavior matrix (,;), as introduced in chapter xiv of Joc. cit. There we 
have considered, among others, the following two extreme cases: complete- 
ly egoistic behavior, characterized by the diagonal matrix 


1 0 0 OMe 


athe a CY CP Ce & 
Ok Oybai 0 eel hearteste tis 
and completely altruistic behavior, characterized by the matrix 
1 1 if Lies eet 
jc te eae 5) 
1 1 REAR ad Eee 1 
An important type of behavior is represented by the matrix 
il 0 0 Ue era 
1 0 0 0 0 
(6) 


© (0 eros? erGe)) 6) les 0 hie Jey ese: eure, el <6 


In this case every individual is concerned with the satisfaction or well- 
being of only one individual. That individual is arbitrarily numbered ‘‘first”’ 
in the matrix (6). This type of behavior is found in patriarchal families, 
where everything is subject to the desires of one man. It is also found in 
autocratic societies in which the autocrat is almost deified and considers his 
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ego superior to everybody else’s. Such societies should not be confused 
with autocratic societies in which the autocrat merely exercises authority 
over the behavior of the rest of the society, but where behavior is directed 
toward a goal different from the mere increase of personal pleasure for 
the autocrat. 

If, as we have considered previously (loc. cit.), the elements k;; of the 
behavior matrix are functions of the social distance s; — s; of the two 
individuals, then both matrices (4) and (5) may be represented as limiting 
cases by a properly chosen expression 


k(s—s’), (7) 


which gives the k’s as functions of s — s’. The general form of those func- 
tions are such that k(x) has a maximum for x = 0, and decreases with in- 
creasing absolute value |x| of x, though they are not necessarily sym- 
metric with respect to x = 0. They thus have the general form of distri- 
bution functions, and it is convenient to introduce the notions of their 
standard deviations. 

If the standard deviation of k(s — s’) is very small, the individual is 
concerned only with himself and with those others who have the same 
social status s, because k(0) = 1. Since the number of individuals is 
finite, while the range of s is continuous, strictly speaking, no two indi- 
viduals will have exactly the same s because the probability of this is in- 
finitesimal. Therefore as the standard deviation o, of k(s — s’) tends to 
zero the behavior matrix tends to the form (4). On the other hand as 
g, — ©, the behavior matrix tends to the form (5); because for ¢, = © 
the individual is equally concerned with all others. 

The elements &;; of the behavior matrix may, however, not only de- 
pend on the social distance s — s’, but also on the geographic distance r. 
Therefore we may more generally consider the k’s as consisting of two 
factors, k, and k,, such that &, is a function of s — s’, while , is a func- 
tion k,(r) of r. Thus 

ROS — Sgr) = he (Cs SR, Cr). (8) 

We shall denote by a, the standard deviation of k,. For very small 
values of o, the individual is concerned only with himself and with his 
immediate “geographical” neighbors, even if individuals located farther 
away have a social status close to his own. For very small values of o, the 
- individual is concerned only with himself, even if there are other indi- 
viduals in close geographic proximity. 

The derivative with respect to s — s’ of the function k,(s — 5’) is 
an even function of s — s’, in the sense that it does not change sign when 
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s and s’ are interchanged though, in general, 
BOCs Ss) 2 hg (Sse (9) 


The function &,(s — s’) decreases when |s — s’| increases. 
The behavior matrix (6) may be considered as a limiting case of a more 


general matrix 


Riu Ri ope, ots: Rin 
Roy Roo tate nea Rone (10) 
kw kyo O00 i Ryn; 


which is characterized by the following relations: 


ki> ka>khg>....> kin (11) 
and 
Rhos — ke eee ane (2) 


That is, all elements of a column are equal, while in each row they de- 
crease from left to right. 

In this general case the maximum concern of each individual is the 
“first”? one, represented by the first column. This first individual also is 
primarily concerned with himself. But each individual has some concern 
about himself, which is, however, less than his concern with anyone to 
the “‘left,”’ but greater than his concern with anyone to the “‘right.’”’ We 
have here a case of complete hierarchy. If we consider the social status as 
decreasing from left to right, then the matrix (10) may be represented for 
a very large number of individuals by making &,(s, s’) a function k,(s’) 
only. Thus, for example, we may put 


RAS) et", (13) 


In this case the concern of an individual with another depends only on 
the social status of that other individual. More generally, we may choose 
k.(s — s’) in such a way that its derivative with respect to s — s’ is an 
odd function of s — s’. For example, we may put: 


k, = ek(s'—s) , (14) 


In this case every individual shows more concern for socially higher in- 
dividuals than for himself and less concern for socially lower individuals. 
But his concern now depends on his own social status, also. A very small 
standard deviation o, of &, in this case reduces the behavior matrix to 
the form (6), while ¢, > © again leads to matrix (5). Thus matrix (5) is 
a limiting case for both cases of even and odd derivatives of k,(s — s’). 
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However, the transition from an even function to an odd one involves a 
discontinuity. The two types of behavior must therefore be considered as 
mutually exclusive different behaviors. 

Both types of behavior can, however, be represented in terms of one 
continuous variable, which we shall call the index of social individualiza- 
tion 7. Let 7 vary in the range from zero to one. For 0 < i < } let the 
society behave according to a matrix characterized by the relations (11) 
and (12), and let o, be a continuous monotone increasing function o,(i) 
of 7, such that o,(0) = 0, o.($) = ©. For } <i <1, let the behavior 
matrix be characterized by k;; > kw and ki; > ki;, and let o, be a con- 
tinuous monotone decreasing function of 7, such that o,(3) = »,¢,(1) = 
0. For « = 0 we now have matrix (6), for i = } we have matrix (5), and 
for z = 1 we have matrix (4). 

The designation “index of social individualization” is justified because 
for z = 0, that is, for matrix (6), the concern of all but one individual is 
for that one. For z = 1 each individual is primarily concerned with him- 
self. For 7 = 4 each individual is equally concerned with himself and 
everyone else. For 7 S 3 the prevailing idea is of the equality of all indi- 
viduals. Even though each individual is concerned primarily about him- 
self, there is no hierarchy present as for 7 < 3. We may, therefore, call 
the quantity E = 27, defined in the range 0 <i < 4, the coefficient of 
equality. 

The desire of an individual to increase, at least to some extent, the 
satisfaction of another depends also on the similarity or difference of the 
satisfaction functions of the two individuals. If the satisfaction functions 
are similar, in other words, if both individuals have similar tastes and 
behave similarly, the altruistic element is likely to be much more pro- 
nounced. For neurobiological reasons discussed in chapter xxii of loc. cit. 
people are apt to dislike behaviors and tastes which are different from their 
own. To take this situation into account we must introduce a third factor 
into k, and write, instead of (8), 


k(s—s’,r,l) =k, (s—s’) k(n) ki), (15) 


where / measures the distance in the functional space between the satis- 
faction function S, of the first individual and Sy of the second. The func- 


tion k,(/) decreases with /. . 
If a1, #2. . . %m are the m possible variables on which the satisfaction 


of any individual in a society may depend, then the satisfaction function 
S of an individual is a function 
MCE ee. (16) 
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of those variables. That function may vary from individual to individual, 
and for some individuals some of the variables may be absent from the 
argument. Assuming a Euclidean metric in the functional space, we may 


put: 
P=f.... f[S. (a, den eee em) 


— S52 (44, Bay oo eee Nd Uae eae cee 


(17) 


where the integration is extended over the whole range of the variables. 

If the standard deviation o of k(l) is very small, then each individual 
will try to maximize only the satisfaction of such other individuals as 
have tastes similar to his own. He will be intolerant of any different tastes. 
If the form of &;(/) is such that for sufficiently large values of /, k,(l) < 0, 
then an individual will be so intolerant of others who have sufficiently 
different tastes that he actually will try to reduce their satisfaction by 
punishment or in some other way. This is a rather common phenomenon. 
When o; = = the individual is equally tolerant of any kind of different 
taste. 

Thus choosing even or odd derivatives of k,(s — s’), as well as the 
values of o,, o,, and o, we can describe quantitatively the types of ethical 
behavior ranging from an adoration of one’s ruler along with complete 
self denial, through an egoistic assertion of the supremacy of one’s own 
person, to complete altruism. Different degrees of tolerance to different 
views, neighborliness, and social “‘snobbishness”’ are all included in this 
general representation. Quite generally we may make k(s — s’, r, l) a 
function of the difference between &,(s — s’, r, 1) and ky(s — s’, 7, 2), in 
a manner similar to the one suggested on page 40 of loc. cit. 

Strictly speaking, the standard deviations o., o,, and o; do not deter- 
mine the situation completely unless the functions k,, k,, and k; are one- 
parametric. The most direct, but also the most difficult, approach to the 
problem is to derive the above functions from some neurobiophysical 
considerations. If the functions thus derived are found to contain more 
than one parameter each, we still can apply the same method of classifi- 
cation. A situation will now be described not only by the a’s, but by a 
larger number of parameters. Pending the development of this suggestion 
we may, as a temporary expedient, postulate some simple one-parametric 
functions. In the following we shall assume that the functions are one- 
parametric. 

Before concluding this section it must be mentioned that other sub- 
divisions of ideas should also be considered in the development of the 
theory. We may classify ideas of behavior according to the content and 
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form of the satisfaction functions. Certain variables may be important in 
one culture and irrelevant in another. This will be reflected in a change of 
the form of S(x1, x2, ... , %m). A variable x; may be considered as impor- 
tant if S(a1, v2, ..., %m) is very sensitive to a change of this variable so 
that |0S/dx;| is very large. A variable «, is unimportant if |AS/dx,| is 
small. A proper classification could be made first in a purely formal way 
and later by ascribing definite meaning to the different variables. 

Another important criterion for classification may be seen in the part 
played by the time element. Individuals may strive to maximize their 
own satisfaction or that of others at the present moment. Or, they may 
be willing to sacrifice now, either for the benefit of future generations or 
for the benefit of a personal “hereafter.” The second case is represented 
by some aspects of Christianity. The first may be found in the contem- 
porary socialistic program of the U.S.S.R. 

We characterize the different degrees of this behavior by an index p, 
which may, for example, indicate the relative frequency of acts intended 
for the benefit of a hypothetical “hereafter” and those intended for the 
gratification of immediate satisfaction. Thus for a mediaeval monk or 
hermit p would be close to unity, while for an epicurean agnostic it would 
be close to zero. 

Last but not least we must consider the classification of scientific ways 
of thinking, a problem that has already been mentioned. The transition 
from the Aristotelian to the non-Aristotelian way of thinking is one case 
at hand. In a previous paper (Rashevsky, 1950) we discussed the “ideal 
scientist” by considering the maximalization of a certain function of the 
knowledge K and of possible information H. It may be possible to con- 
struct such a one-parametric function of K and H that for different values 
of the parameter the maximalization leads to a greater or lesser preference 
for either K or H. In this way a continuous transition from Aristotelian to 
non-Aristotelian thinking may be quantitatively described. 

Il. Effects of Ideologies on Social Structure. In loc. cit. we saw that the 
distribution function f(s) of the social status s is determined by two other 
functions. One is the probability that an individual of ability a and social 
status s will change his social status by an amount A when he interacts 
with another individual of ability a’ and social status s’. This function 


is designated by 
DCG Ee 8,5 4 Adi (18) 


The other function determines the relative frequency of interaction of 
two individuals in terms of their respective social statuses. In Joc. cit. it 
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was designated by K(s, s’). Here we shall denote it by 
RSs ys Ge) 


It is readily seen that both functions depend on, and perhaps are com- 
pletely determined by, the functions k.(s — s’), k-(r), and k,(l). The function 
R(s, s’) is in many cases of the form R(s — s’). In that case the smaller 
its standard deviation cz, the more often will interactions occur only be- 
tween individuals of nearly equal social status. Since the interaction of 
two individuals usually implies some degree of mutual concern, there- 
fore, in general, o will increase with o, and o;. We should also expect the 
standard deviation o, of p(a, a’, s, s’, A) with respect to A to decrease 
with decreasing os, ¢,, and o;. A poor or socially low individual interacting 
with a rich or socially prominent person, who is concerned only with him- 
self, is not likely to improve his social or economic position as a result of 
such interaction. Also, if 7 is great and o, small, the interaction may be 
only an indirect one—say, by mail—and is likely to be less effective. 

As before, we may again either attempt to derive the functions p(A) 
and R(s — s’) from the functions k,, k,, and ki; or, as a temporary ex- 
pedient, we may postulate one-parametric functions for p(A) and R(s — s’) 
and postulate simple relations, e.g., linear ones, between o, and cp on 
one hand, and oa,, o,, and a on the other. Since the vanishing of either 
os, 0,, or o; implies vanishing of o, and og, therefore the latter two may 
be linear in each of the former variables, but the three former should 
enter in the form of a product. Thus we may put 


Tr =CROs0;01 ’ (20) 
Op = CaO 5or Cr (21) 
where C and C, are constants. 
To include the possibility of pure chance we may put, generally, 


or =CROo,0;+ 0%, (No, ni) 5 (22) 
op =Cyosa,a1+ 0% (No, ni) 5 (23) 


where of and o, depend on the total population No and on other param- 
eters 7;, such as area occupied by the population, means of communica- 
tions, etc. 

If we can solve the functional equation (21) of chapter viii of loc. cit. 
then we can express f(s) in terms of the parameters o., ¢,, and o; which 
determine the classification of ideas. 


Again it may be possible to roughly classify all functions f(s) by their 
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standard deviations. The problem therefore arises as to whether it is 
possible to express approximately the standard deviation of f(s) in terms 
of the standard deviations og and o, without actually determining f(s) 
completely. 

III. Fundamental Equation of the Dynamics of Ideologies. Let us now 
apply the reasoning of chapter xxii of our book (Rashevsky, 1951a) to the 
case of 2 mutually exclusive behaviors. Following a suggestion of H. D. 
Landahl’s (1941) which we have developed somewhat further in a recent 
paper (Rashevsky, 1951c), we now shall consider the m quantities 


ge (hx). (24) 


The quantity ¢; now determines the probability that behavior R;, will 
be exhibited. If g; = 0, then R; is as probable as any one of the other 
n — 1 behaviors. It is then equal to 1/n. As 9; > &, the probability 
P; of behavior R; tends to 1. If 9; < 0, and tends to — ~, then P;— 0. 
The actual expressions for P; have been given elsewhere (Rashevsky, 
1951c). 

If 1 is very large we may substitute integration for summation. Let 7 
be a continuously varying parameter, each value of which corresponds 
to a particular choice of behavior. For example, in some cases » may 
be identified with the index of social individualization i. Or 7 may 
stand for o, or o;. Let 7 vary in the interval (m, n2). For every value of 
there is a corresponding value e(7). Denoting by 7, the value of 7 for which 
v(n) has a highest value, we now put 


e (nm) =€(nm) —€é(n) (2:5) 
where 
1 19 
é = 26 
Oe GONG (26) 


and the integration is taken over all values of 7 except 1m. 

If the total ideology or behavior is characterized by several parameters 
ni, then the e’s and g’s are functions of all of these parameters. Without 
loss of generality we shall, for simplicity, confine ourselves here to the: 
one-dimensional case. i 

With the definitions given above we now may apply the same method 
as for the choice between two different ideas or behaviors. The prob- 
ability P(n) of an ideology 7 may now be given by 


eh ere ere (27) 
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where o is a constant, or by: 


reg AES, (28) 


K 


where x, is as defined in Rashevsky (1951c). 

It may at first glance seem more appropriate to consider different 
values of n as corresponding to different intensities of the same behavior. 
However, considering the actual inability of most persons to see essential 
similarities between only slightly different situations and to realize that 
the differences are of quantitative rather than qualitative nature, it seems 
justifiable to use the approach suggested above. A good example of this 
is the bitter antagonism shown each other by rival political parties whose 
political ideas actually differ but slightly. 

The approach suggested above is justified moreover from a neurobio- 
physical point of view. The mathematical biophysics of discrimination of 
intensities (Rashevsky, 1948, chap. xxiii) leads to the concept that quali- 
tatively similar stimuli of different intensities result in the excitation of 
spatially different regions in the central nervous system. Thus even sev- 
eral different intensities of the same type of behavior lead neurobio- 
physically to the same model as do several spatially different mutually 
exclusive stimuli. 

The neurobiophysical picture will, moreover, help us to establish the 
equations which govern the shifts from one value of 7 to another, equa- 
tions which were very tentatively suggested in our book (Rashevsky, 
1951a). 

Consider a structure consisting of 7 parallel, mutually cross-inhibiting 
pathways (Rashevsky, 1948, chap. xxxi). Let the ith pathway be excited 
at the afferent end by an amount e; of the excitatory factor. Let the 
parameters of the system be such that when all e;’s are equal, all the re- 
sponses at the efferent ends are just inhibited. If a single e, exceeds the 
average value e of the other e,’s, then the reaction R; is produced. For 
very large this picture leads to the equations (25)—(28), if each reaction 
_R; represents the behavior corresponding to a particular value 7; in our 
quasi-continuous scale. 

The excitations e,’s are considered to be of a central origin, due to some 
endogenous factors, just like ¢, and & in the case of only two alternative 
behaviors. However, let the centers which produce the e’s be connected 
by cross-excitatory pathways along which the intensity of excitation de- 
creases with distance, as was studied previously (Rashevsky, 1948, pp. 
381-82; Landahl, 1938). Thus if a given e; is increased, then the e’s in 
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the immediately adjacent centers will also increase, this increase diminish- 
ing with distance from the given center. For definiteness let us consider 
the case in which the decrease of the cross-excitation with distance is 
exponential. If we visualize all the m pathways arranged uniformly in a 
linear array, then the actual distance measured along this array between 
two points will be proportional to the difference 7; — 7, ,where n: and 
are the values of » which correspond to those two points respectively. 

From the above it follows that if all e’s remain equal to each other, 
except one, €(7m), which for any reason increases, then the values of e(n) 
in the neighborhood will also increase, so that e(7) will now be repre- 
sented by 


Gee Gale fe (29) 


where €p is the level of excitation of the other centers and wy? is a constant. 
If uw? is sufficiently large, the distribution of e(n) will represent a sharp 
peak in the neighborhood of 7. 

From equations (25), (26), and (27), together with equation (29), it 
follows that the probability P(7) is highest for 7 = nm». Without imitation 
the largest percentage of all individuals will exhibit it. In the presence of 
imitation an absolute majority will exhibit behavior 7,,. The actual per- 
centage of that majority is obtained from the previously established equa- 
tions (Rashevsky, 1951a), in which we substitute ¢(7,,) for ¢. 

As we have seen, the choice of a particular behavior 7; results in a par- 
ticular social distribution function f;(s). The latter in turn determines 
the rate at which the chosen behavior will be either enforced or weakened 
by learning. Whether it is reinforced or weakened depends (Rashevsky, 
1951a, chap. xxii) on the sign of the expression 


f KOM) as, (30) 
0 


where X(s) is positive for those values of s for which the behavior results 
in satisfaction, and negative for those values for which the behavior is 
unsatisfactory. We may, for example, put: 

N(s) = 225 —(? , (31) 


where a and @ are constants. If A(s) is to be positive for s = 1, then 
a? > 6, If a? = 267, then 

N(s) =a? (s— 3) (32) 
so that d(s) = 0 for s = 4. It seems more plausible that d(s) = 0 for 
s <4, which implies a? > 6’. In line with the discussions of the 
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satisfaction function (Rashevsky, 1951a) it may be more logical to put 
h (s)*= log a?'s (33) 
with 

oS, (34) 
Generalizing the concepts developed in chapter xxii of our book (Ra- 
shevsky, 1951a) we may consider that a negative total satisfaction reduces 
the value of ¢(m) by reducing €(7m). Due to the cross-excitatory connec- 

tions, this will also reduce the adjacent e(»)’s. 
Without loss of generality we may put 7, = 0 and let » vary in the 
interval (0, 1). Let the society initially have a preference for behavior 
and ideology 7 = 0, and let the value of e(7,,) at the time ¢ = 0 of adop- 


tion of that ideology be 
e(nm) =&. (35) 


If that ideology or behavior results in an average dissatisfaction, then é 
will decrease with time, and at any time ¢ > 0 we shall have 
e(n) =@q+é& (¢) em. (36) 
Eventually é(£) will become negative and e(n) will have its smallest 
value at » = 0 and largest at » = 1. That largest value is equal to 
éo + & () ew me. (37) 


The value of g(1) now increases. At any moment fit is given, according 
to equations (25), (26), and (37) by 


y (1) =¢£ (2) et — t (1) ff ewndy =— E w) Shas (t) ae enn. (38) 
0 M M 


If, as we shall assume, yw? > 1, then, approximately, 
t 
e(t) = EP so. (39) 


As soon as ¢(1) > ¢*, where ¢* is the threshold determined by the 
slope of the F(y) curve and the slope of the ay line (Rashevsky, 1951a, 
chap. xii), society will switch into ideology and behavior 7 = 1. This will 
happen when 


—&() =n%¢". (40) 


If this behavior results in a net dissatisfaction, then the value of e(1), 
which at the moment of the switch was € + ge” ~ €9, will begin to 
decrease, and at every moment ¢ be given by 


ato, ¢) <0, (41) 
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if we now put as the new origin of time ¢ the moment when ideology n = 1 
is adopted. 

Because of the excitatory cross connections this will result in an addi- 
tion to every e() of the (negative) amount 


¢ (1) enn), (42) 


At the same time the original inhibition, which reduces €, is now de- 
creasing, due to “forgetting” (Rashevsky, 1951a, p. 180). This decrease 
may be assumed to proceed exponentially with time. Therefore, we now 


have 
e(y) =@ — Mg e vte= Hn — C(t) ew —n) , (43) 


The function ¢(7) has a maximum at 
1 
max = 7 3 lu? — yf — log 6 (t) +log p2¢*. (44) 


This maximum is somewhere between 7 = 0 and 7 = 1, and shifts 
with time. It is readily seen that it also becomes more pronounced. By in- 
troducing expression (44) into (43), and subtracting from this the inte- 
gral of (43) with respect to 7 from 0 to 1, we find the expression for 
¢(nmax) aS a function of ¢. Setting this equal to ¢* we find the moment #, at 
which the society will switch to a new behavior. By introducing this value 
of t for ¢ into expression (44) we find the max which describes the new 
ideology and behavior. From there on we again repeat the procedure, 
remembering that now expression (42) is also to be multiplied by e~”. 
In addition, since the new behavior lies between 7 = 0 and 7 = 1, we 
shall have to study both branches of the curve (29). 

We thus have given a possible analytical treatment to the problem out- 
lined on pp.. 180-81 of our book (Rashevsky, 1951a). 

There is, however, another possible way for society to shift from one 
ideology or behavior to another. If expression (30) is positive, which is 
equivalent to a positive average satisfaction, then the value of ¢(n) for 
the chosen behavior y will increase. In the absence of other factors this 
would make the behavior 7 permanent. The fraction of individuals ac- 
cepting it would tend, with time, to one. 

We must remember, however, that the quantities z, ¢,, and o:, which 
may stand for 7, vary in general with increasing intellectual development 
of society. They will depend on the knowledge K, and especially on the 
amount of communication transmitted from individual to individual in a 
society. The greater the density of population, the more likely is an in- 
dividual to come in contact with others and, therefore, the more likely he 
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is to know other individuals well. This better knowledge is likely to in- 
crease o;, o,, and o;. Thus, even in the absence of any satisfaction or dis- 
satisfaction, that is, even if expression (30) is zero, a slow continuous shift 
of 7 will occur. Mathematically this can be taken into account by making 
the quantities which correspond to 7 in equation (29) be slowly varying 
functions of time at every step of the previous discussion. 

Thus when a particular ideology or behavior 7 is chosen, two things 
may happen: 

If the chosen behavior results in discontent, then it will eventually be 
abandoned for one which is remote on the scale of behavior. The transi- 
tion to the new behavior or ideology will be sudden when the threshold 
y* is exceeded. 

If the chosen behavior is, on the average, emotionally neutral, or even 
pleasant, then a slow shift in the chosen 7 will occur. From a behavior 
characterized by a value of 7 a change will be made to one characterized 
by mm + An, where Ay is a small increment. The change will be, strictly 
speaking, discontinuous, in the sense that the majority of individuals who 
exhibited behavior 7 will suddenly begin to exhibit behavior 7m + An. 
Because of the smallness of An, however, the total process will appear to 
be continuous. 

From equations (30)—(34) we see that the total dissatisfaction in- 
creases with increasing steepness of the distribution function f(s). If, as 
a measure of that steepness, we take the standard deviation o; of f(s), 
then we see that for small values of o; sudden changes from one behavior 
to another will occur. For sufficiently large values of o; the transitions will 
be continuous. This agrees with the results at which we arrived previously 
(Rashevsky, 1951a, chap. xxvi) by a different method. 

The previously developed equations which govern imitative behavior, 
together with equations (25) and (26), give us the percentage X of the in- 
dividuals who exhibit the chosen behavior (Rashevsky, 1951d). The per- 
centage VY now includes all the individuals who exhibit all the non-7,, be- 
haviors. It is interesting to inquire how the other behaviors will be dis- 
tributed within the fraction Y. It is natural to assume (and an attempt to 
prove it may be worthwhile) that the relative portion of individuals among 
the fraction Y who exhibit different non-7, behaviors is proportional to 
the corresponding values of ¢(7) for 7 ~ mm. Or it may be that of all the 
behaviors 7 ¥ 7, the one is preferred which has the highest value of (7) 
and the preference is determined by a mechanism similar to the preference 
for mm among all 7’s. 

Before any applications of the theory outlined can be made it will have 
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to be developed considerably. By way of an illustration as to what eventu- 
ally may be expected from such a theory we may very tentatively suggest 
the following interpretation of some historical processes. 

Toward the end of the Roman Empire behavior was characterized by 
a matrix rather similar to (6). The index 7 of individualization was near 
zero for the majority of the population. The subjugated people and the 
vast number of slaves worked for their masters, whose number was rela- 
tively small. 

Due to factors discussed in section I of this paper the ideology of mono- 
theism then spread in the form of Christianity. If we classify the latter in 
our system of ideologies, we see that with its idea of the brotherhood of 
man it corresponds essentially to matrix (5), with = 4ando, = 0, = ©. 
Another characteristic of Christianity was the development of asceticism 
with belief in reward in the hereafter; hence with a value of p near unity. 
Of the above characteristics only the second prevailed in practice. Soon 
after the beginning of Christianity the idea of universal brotherhood be- 
came non-existent in practice. The large value of p is likely to have been 
caused by the highly unsatisfactory position of large oppressed masses. 
Those masses could relieve their dissatisfaction only in one of two ways: 
by gaining actual control of the society and improving their situation, or 
by reducing their dissatisfaction through reducing their desires and 
through hopes for a better hereafter. Physically the two ways lead to very 
different results, but psychologically the results are similar [cf. P. Soro- 
kin’s definition of freedom (1937)]. 

Thus at the beginning of the Middle Ages ideology and behavior were 
characterized by 7 = 0, 0, X 1, 0; « 1, p ~ 1. This, as we have seen, 
is likely to lead to a very steep f(s), that is, to a very small o;. Though no 
quantitative data on f(s) in the Middle Ages are available, the social 
stratification is known to have been very pronounced. This bears out our 
conclusion. 

But a small value of o; should have resulted in an eventual sudden 
transition to an ideology with z = 1, p ~ 0, 0, ~ 1, 0. ~ 1. This sudden 
transition should, however, have been preceded by a relatively long period 
with an apparent lack of change. As we have seen before (Rashevsky, 
1951a) the gradual shift in ¢, or, in our case, in ¢(7) results in a parallel 
displacement of the curve F(y) (Fig. 1) to the right. If the slope of the 
straight line is very small, or if the initial ¢g is sufficiently positive so that 
the intersection with the F(y) curve is very near the asymptote, then a 
displacement of that curve will not materially alter the value of X, until 
the curve reaches a position somewhere in the neighborhood of the broken 
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line. From there on X begins to decrease appreciably, until the position 
indicated by the dotted line is reached. At that point X changes discon- 
tinuously to a value which is near zero. 

The period between the position of the full line and that of the broken 
line may be considered to correspond to the Middle Ages proper. The 
position of the broken line may be considered as the beginning of the 
Renaissance, while the position of the dotted line may be considered to 
correspond to the quasi-discontinuous revolutionary changes of the end 
of the 18th Century and the beginning of the 19th. Those changes, as 
we have seen, should have resulted in ideologies which correspond to 


Ficure 1 


p<<«<1,72~1. But ¢ = 1 corresponds to matrix (4), which is charac- 
teristic of extreme individualism. Economically matrix (4) corresponds 
to what is usually designated as “free enterprise” (Rashevsky, 1951a, pp. 
146-7). A small p means, in common parlance, a more materialistic out- 
look on life, with the idea of any “hereafter” taking a very subordinate 
place. But such changes have been characteristic of the late 18th and 
early 19th centuries. 

While the values of 7, ¢,, and o, affect the form of f(s), they do not de- 
termine f(s) completely. A number of other factors also affect f(s), 
especially if s denotes the economic status. Therefore after the idea of 
individualism established itself in different countries, its further fate 
varied from country to country. In western Europe with a high density 
of population and high industrialization o; was relatively large. Hence 
further shifts in 7 and in the other parameters took place quasi-continu- 
ously. The results of previous studies (Rashevsky, 1951a, chapters xiv, xv) 
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indicate that a minimum total dissatisfaction is likely to occur in a society 
characterized by matrix (5), that is, by i = 4. Hence we should expect a 
quasi-continuous shift toward such a behavior in western European coun- 
tries. On the contrary, in a country with little industrialization and a high 
degree of total dissatisfaction, like czarist Russia, a sudden transition 
to extreme socialism was to be expected. 

Some of the conclusions outlined here may in principle be checked by 
gathering qualitative or semi-qualitative data on socio-historical develop- 
ments of the type, for example, compiled by P. Sorokin in his Social and 
Cultural Dynamics (1937). Among other things Sorokin attempts to fol- 
low by centuries the fluctuation of such qualities as the relative frequency 
of appearance of religious and secular subjects in art, or the relative fre- 
quency of description of different social classes in art. Such figures, in- 
exact as they are, may be indications of the relative values of X and Y, 
where X may stand, for example, for the fraction of individuals who 
exhibit characteristically mediaeval behavior, with o, «1, ™ <1, 
p ~ 1; while Y represents the fraction of individuals who exhibit the op- 
posite behavior which is characteristic of modern times. Before making 
use of such data, even if they were accurate, an important question must 
be decided. Do the values of X and Y as given by such data represent the 
actual values for the society at a given period, as determined by mass 
behavior and given by ¢ + y, or do they reflect the innate values of ¢, 
which would manifest themselves in the absence of imitation? At first 
glance the former possibility seems to be the natural one. However it is 
possible that the artist reflects the inner tendencies of his fellow men more 
than their overt behavior. As we have seen (Rashevsky, 1951a), the aver- 
age @ may have already shifted well toward the opposite behavior, yet the 
majority of society may still exhibit the given behavior. 

In order merely to illustrate what could be done with such data, let us 
assume the second situation. In that case X is determined not by the in- 
tersection of the curve F'(W) with the straight line ay, but by the ordinate 
of the curve F(¢), where ¢ is the average value of ¢ at a given moment. 
We have 


X=7ll+F le) | (45): 


and, approximately, 


F (¢) =G (2), (46) 


because the quantities x», (Rashevsky, 1951c) are very small for large 
and may be neglected. 
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Let us assume that ¢ varies linearly with time: 
G=v(t—b), (47) 
and that positive ¢ corresponds to “mediaeval” behavior. We now intro- 


duce the new variable 


gy =2= (t-h). (48) 
Then, putting 
Ki = (49) 
we have, from (46) and (47): 
PG=h) 2G (—). (50) 


From the values of X at various times, we find 
F(t—t) =2X—1 (51) 


and then we may calculate x’ and éo. In Figure 2 we have plotted the values 
of F computed from Sorokin’s data on the relative percentage of religious 
and secular art for central Europe (Sorokin, 1937, Vol. I, p. 381), marked 
by squares; data on same frequencies for Europe as a whole (ibid., p. 
382), marked by asterisks; and data on relative frequencies of rep- 
resentation of different social classes (ibid., p. 486), marked by tri- 
angles. The full line represents the G-curve, which fits those data as well 
as is possible in view of their inaccuracy. Sorokin’s data (Fig. 2) gives us, 
if we use 100 years as the unit of time, x’ = 0.5, fp) ~ 1650 a.p. We see 
that o’ = 0 and hence ¢ = 0 at about 1650 a.p. Since the actual revolu- 
tionary transition toward negative ¢’s occurred at the end of the Eight- 
eenth Century (French Revolution), therefore this supports our assump- 
tion that Sorokin’s data gives us @ rather than g@ + y*. For ¢ > fy the 
F curve, as seen from (50), is merely displaced to the right. At approxi- 
mately 1800 a.D. it must have been in the position shown by the broken 
line of Figure 2. At that time the sudden transition occurred; hence it was 
then that the straight line ay became tangent to the curve. Since the 
transformation (48) does not alter the ordinate of the point of tangency, 
but merely changes the units of g so that numerically ¢ becomes equal to 
t — to, by drawing a tangent to the F curve in Figure 2 from the point 
t — to = 0, we find the value of F and hence of X at the moment of the 
revolution. The value is approximately X = 90%, which is quite plausible. 
After the revolution we have X ~ 0.02, Y ~ 0.98, which is again plau- 
sible. In principle those values could be determined by appropriate quanti- 
tative historical research. 
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One more illustration may be given. During the Middle Ages X was of 
the order of 98%, Y of the order of 2%. As we have seen elsewhere (Ra- 
shevsky, 1951b), the scientific and technological productivity of a society 
should be approximately proportional to Y. If this is the case we should 
expect approximately a 50-fold increase in that productivity. Since the 
French Revolution P. Sorokin (1937) finds the approximate number of 
discoveries and inventions made in western Europe between 600 a.p. and 
1800 a.D. to be 3,000. From 1800 a.p. to 1900 .p. he finds that number to 
be approximately 8,500. Per unit time this is an increase of 34 times. 

We again must emphasize that all the above is offered only as an illus- 
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tration of how mathematical reasoning could be applied to historical proc- 
esses. Before this actually can be done a great deal of further development 
of the theory will have to be made. Moreover only if such a development 
stimulates quantitative historical research, which now is practically non- 
existent, can any applications be made. 

Important quantitative questions arise in connection with comparative 
studies of different cultures. Why did the West develop earlier industrially 
than the Orient? What factors determine the absolute rates of develop- 
ment? Those questions we do not attempt to discuss here, but hope to do 
so in future publications. 

The author is indebted to Dr. George Karreman for a critical discussion 
of the paper. 
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AWARDS FOR POST-DOCTORAL STUDY IN STATISTICS 
AT THE UNIVERSITY OF CHICAGO 


The Committee on Statistics (a department) of The University of 
Chicago has established, under a five-year grant from the Rockefeller 
Foundation, a program of Post-doctoral Awards to provide training and 
experience in statistics for scholars whose main interests lie outside that 
field. There will be three Awards per year, to holders of the doctorate or 
equivalent in the biological, the physical, and the social sciences. Each 
Award will be $4000 or slightly more, office space will be provided, and 
$600 to $1000 will be available for clerical, computational, and research 
assistance. There will be no tuition charges. 

The purpose of the Awards is to give statistical training to a few scien- 
tists who may be expected to employ it both to the direct advance of their 
specialties and to the enlightenment of their colleagues and students by 
example, by consultation, and by formal instruction. The development of 
the field of statistics has been so rapid that problems of communication 
are a serious obstacle to its full exploitation. The amount and quality of 
instruction available to young students is constantly increasing, but there 
is a real need, which these Awards seek to fill, for making appropriate 
instruction available to already established scientific workers who give 
promise of immediate applications of statistics to their special fields. 

Recipients of the Awards must have received the doctor’s degree prior 
to commencing the program, except in the case of recognized research 
workers whose experience and accomplishments are clearly the equivalent. 
Candidates whose mathematical preparation includes less than the usual 
sophomore year of calculus, or its equivalent, will not ordinarily be con- 
sidered, but previous training in statistics is not required or expected. 
Candidates having under way research programs in their own fields will 
be preferred, and the department of The University of Chicago concerned 
with a candidate’s specialty will be asked to participate in evaluating his 
application. Recipients must spend eleven months studying statistics at 
The University of Chicago, and will be expected to pursue a number of 
regular courses. 

Applications, or requests for further information, should be sent to: 
Committee on Statistics, The University of Chicago, Chicago 37. Applica- 
tions for the academic year 1952-53 should arrive by April 1, 1952. 
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